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Abstract 



The top-pair ti production cross section close to threshold in e~^e~ collisions is strongly affected 
by the small lifetime of the top quark. Since the cross section is defined through final states con- 
taining the top decay products, a consistent definition of the cross section depends on prescriptions 
how these final states are accounted for the cross section. Experimentally, these prescriptions 
are implemented for example through cuts on kinematic quantities such as the reconstructed top 
quark invariant masses. As long as these cuts do not reject final states that can arise from the 
decay of a top and an anti-top quark with a small off-shellness compatible with the nonrelativis- 
tic power-counting, they can be implemented through imaginary phase space matching conditions 
in NRQCD. The prescription-dependent cross section can then be determined from the optical 
theorem using the e~^e~ forward scattering amplitude. We compute the phase space matching 
conditions associated to cuts on the top and anti-top invariant masses at next-to-next-to-leading 
logarithmic (NNLL) order and partially at next-to-next-to-next-to-leading logarithmic (N^LL) or- 
der in the nonrelativistic expansion and, together with finite lifetime and electroweak effects known 
from previous work, analyze their numerical impact on the tt cross section. We show that the phase 
space matching contributions are essential to make reliable NRQCD predictions, particularly for 
energies below the peak region, where the cross section is small. We find that irreducible back- 
ground contributions associated to final states that do not come from top decays are strongly 
suppressed and can be neglected for the theoretical predictions. 
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I. INTRODUCTION 



The measurement of the hne shape of the total cross section a{e~^e — tt) for top quark 
pair production at energies around the top-antitop threshold {\/g^ ~ 350 GeV) constitutes a 
major part of the top physics program at a future linear collider (LC). The rise and the form 
of the cross section allow for precise measurements of the top quark mass rrit in a threshold 
mass scheme [l|, the top quark width Ff, the top- Yukawa coupling i/t and the strong coupling 
as- In view of the expected experimental precision at the LC [2|, theoretical uncertainties 
for the predictions at the level of da/cr ~ 2 — 3% in the peak and the continuum region are 
desired In the energy region below the peak where the cross section is becoming tiny 
the theory error should not exceed the level of around 5 fb 0]. 

From the theoretical perspective the QCD dynamics of the top quark pair for threshold 
energies is quite nontrivial since the small relative velocity of the top-antitop pair v 1 
leads to a proliferation of physical scales that need to be accounted for: the top mass 
rrit, the relative three-momentum p ~ rritV and the nonrelativistic kinetic energy E ~ 
mtv"^ of the top quark pair. In the standard QCD perturbative expansion singular terms 
oc (as/v)"" and oc (oslnw)" arise from the ratios of these scales. Since the Coulomb-like 
dynamics enforces the power-counting v ~ asirritv) the proper treatment of these terms 
requires resummations within a power-counting framework with a simultaneous expansion 
in w ~ as <^ 1. This can be achieved with nonrelativistic QCD (NRQCD), a low-energy 
effective field theory (EFT) of QCD that separates the different quantum fluctuations that 
are relevant for the kinematic situation of heavy nonrelativistic quark pairs. Within the 
fixed-order approach, which achieves a systematic summation of terms cx: complete 
NNLO predictions (i.e. n+m < 3) have been made For available NNNLO results we refer 
to Refs. |3-[l3 • Using a proper low-scale short-distance threshold mass scheme with cut-off 
scale R ~ rritas fl3l-16| (which avoids the pole mass renormalon) the energy where the cross 
section rises is stable in perturbation theory. However, the fixed-order predictions suffer 
from large normalization uncertainties at the level of 10-20 % which indicate potentially 
large logarithmic terms. These normalization uncertainties are particularly problematic for 
measurements of the top width and the Yukawa coupling yt. This problem of potentially 
large logarithmic terms is addressed in renormalization group improved NRQCD calculations 
of the cross section, which account for a systematic summation of terms oc a"?;™ln^w. 
Renormalization group improved QCD results are fully known at NLL order (i.e. n + m — 
k < 2). At NNLL order (i.e. n + m — k < 3) all ingredients are known except for the 
NNLL renormalization group evolution of the Wilson coefficient of the leading order top 
pair current Partial results for the NNLL order anomalous dimension of the current 

have been computed in Ref . 0, [20[ . At the present stage renormalization group improved 



results have a normalization uncertainty of 6-10 % [19|, |21 . 

Electroweak effects and in particular the top quark decay play an equally important role. 
Already at leading order it is important to account for the top decay width since it widens 
the top-antitop bound state resonances and turns the threshold cross section into a smooth 
lineshape. Despite this fact the determination of subleading electroweak and finite lifetime 
effects have received somewhat less attention in the literature in the past. Theoretically, for 
predictions of an inclusive cross section the top width acts as an infrared cutoff for the top 
energy and thus allows for a perturbative computation for all threshold energies. We have 
Ft ^ {GF/8V2n)mf ^ 1.5 GeV, which scales like g'^rrit ~ gjrrit where gi and g2 are the U(l) 
and SU(2) couplings, and we also find numerically that Tt ~ mto;^. With (^i ~ (?! ~ Q^s it is 
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therefore natural to adopt the power-counting 



v^as^ aj^l < 1 (1) 

when making perturbative theoretical predictions. 

Electroweak interactions are responsible for various effects that can be categorized into 
four classes: (a) "Hard" electroweak effects, which includes hard, point-like corrections 
related to the tt production mechanism by virtual photons and Z exchange, as well as 
electroweak corrections to the hard matching conditions of the NRQCD operators and po- 
tentials, (b) electromagnetic effects for the luminosity spectrum of the e~^e~ initial state, 
(c) electromagnetic corrections to the low-energy nonrelativistic dynamics of the tt pair and 
its decay products and (d) effects related to the finite top quark lifetime. The corrections 
from class (a) can be determined by standard methods through matching for top quarks 
in the on-shell limit and are real numbers. Up to NNLL order they have been discussed 



and implemented in Ref. [22[. Earlier work can be found in Refs. [23|, |2J]. The QED beam 
effects from class (b) are taken into account by a convolution of the "partonic" cross section 
with an on-shell e+e~ pair in the initial state with the collider's luminosity spectrum. The 
luminosity spectrum accounts for initial state radiation, the accelerator-dependent beam 
energy spread and the beam-strahlung, and effects coming from beam-beam interactions. 
Since the luminosity spectrum is for the most part determined either experimentally or from 
experimental simulations we do not consider it any further in this work. The photon in- 
teractions of the ti pair belonging to class (c) are quite similar to the gluonic corrections 
and can be incorporated in the same way into NRQCD through potentials and ultrasoft 
interactions. The most important of these effects are the QED corrections to the Coulomb 
potential contributing at NLL order. 

The finite lifetime effects of class (d) are the main purpose of this work. A proper treat- 
ment of instability effects entails that for the definition of the cross section, final states 
compatible with the top decay chains are accounted for. Since this necessarily also includes 
final states that do not arise from top decays, but are experimentally indistinguishable, the 
cross section - and potentially also the theoretical methods to compute it - are dependent 
on what selection prescriptions are employed. In general, as long as selection cuts do not 
reject final states that can arise from the decay of a top and an anti-top quark with a 
small off-shellness compatible with the power counting of a nonrelativistic (anti-)top quark 
propagator^, they can be implemented through imaginary contributions to the Wilson co- 
efficients of the NRQCD operators. This means that the top decay products are integrated 
out together with the information on the selection cuts. The resulting NRQCD Lagrangian 
is formally non- Hermit ian. The cross section can then be obtained via the optical theorem. 



i.e. from the imaginary part of the e~^e~ — )■ e~^e~ forward scattering amplitude [25|-|27|. In 
this work we adopt this inclusive approach. 

One can distinguish two different types of imaginary contributions to the NRQCD Wilson 
coefficients: (1) contributions arising from cuts of full theory diagrams through top decay 
final states and (2) contributions describing the selection prescriptions, which we will also 
frequently refer to as cuts in the following. ^ The type-1 contributions describe decays of 
(anti)top quark modes that are propagating in NRQCD. In the matching procedure they 



^ Wc call such selection prescriptions "inclusive" throughout our work. 

^ The different meanings of the word "cut" used frequently in this work should be clear from the context. 
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arise from cuts of full theory diagrams through (anti)top decay final states. The leading 
type-1 NRQCD term is the well known on-shell width contribution, 



5£ = 5^^t^r,^p + 5^xi^r,xp, (2) 
p p 

which in the matching procedure comes from the bW cut of the full theory top quark on-shell 
self-energy, and which contributes at leading-log (LL) order in NRQCD. Here Fj is the top 
quark on-shell width, and ipp and Xp represent Pauli spinor field operators that destroy top 
and antitop quarks, respectively. This term leads to the NRQCD (anti)top propagator of 
the form 

(3) 



-n — P -4- -i^t 



and can be readily implemented into computations for stable top quarks supplemented by 
the replacement rule E ^ E + iTt, where E = ^/s — 2mt is the cm. energy with respect 



to the two-particle threshold [25|, |26|. The type-1 contributions to the Wilson coefficients 



up to NNLL order and neglecting the width of the W-boson were determined in Ref. [27 
As a new higher order feature these include interference contributions from double-resonant 
(e+e" tt ^ hhW^W~) and single-resonant (e+e~ tbW~ ,tbW^ -)■ bbW~^W~) am- 
plitudes as illustrated in Fig. [H Electroweak gauge-invariance is maintained during the 
matching procedure by the inclusion of the imaginary contribution of the top quark wave 
function renormalization Z-factor arising from the bW intermediate state in the top self- 
energy of the full theory. 

An important new theoretical aspect is the emergence of imaginary anomalous dimen- 



sions caused by UV-divergences in the NRQCD tt phase space integrations [27[. These 
UV-divergences originate from the Breit-Wigner behavior of the top and antitop propaga- 
tors of Eq. ([3]) which, upon being cut, lift the on-shell dispersion relation pq = p^/2mt. 
In the forward scattering amplitude these propagators allow (anti)top intermediate states 
contributing to the NRQCD cross sections which have arbitrarily large invariant mass. 
Once f ^-suppressed NNLL order operators are inserted, this behavior causes UV-divergences 
which are compensated by imaginary counter-terms associated with {e~^e~){e~^e~) forward- 
scattering operators. The resulting NLL order imaginary anomalous dimensions of the 
Wilson coefficients of the {e~^e~){e~^e~) forward-scattering operators sum \nv terms from 



the top decay phase space. It was shown in Ref. [27[ that the type-1 contributions to the 
imaginary parts of the NRQCD Wilson coefficients are numerically important for the nor- 
malization of the cross section and for the correct prediction of the cm. energy of the tt 
quasi-resonance peak. The existence of phase space divergences in NRQCD inclusive cross 
section computations for top pair production at threshold indicates that - due to the finite 
lifetime - predictions beyond the leading order approximation need additional short-distance 
information to be defined unambiguously. This short-distance information is provided by 
the experimental selection criteria for the final states that are accounted for the cross section 
determination, and also incorporate background contributions related to diagrams without 
an intermediate ti pair.^ In NRQCD this phase-space short- distance information is incor- 
porated in the imaginary contributions to the Wilson coefficients of type-2 we have already 



Also outside the framework of effective tlieories such selection criteria are in general necessary, and 
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FIG. 1: Interference of the double-resonant diagram with a single-resonant one, as described in 
the full and effective theories. The (grey) blob in the tt production vertex in the EFT diagram 
represents the NNLL matching condition which results from integrating out the lower bW~ loop. 
The cut through this vertex represents the contribution of the absorptive part of the matching 
condition in the optical theorem. The (red) dashed line means we extract the imaginary part of 
the forward scattering amplitude or, equivalently, that we perform the phase space integration over 
the particles in the cut. Note that we use double lines for representing the top quarks. 



mentioned above. We call the computation of these contributions the "phase space match- 
ing" procedure. Many different feasible inclusive cross section definitions can be devised. 

In this work we determine and analyze the type-2 imaginary contributions to the NRQCD 
Wilson coefficients at NNLL and N^LL order for cuts on the invariant masses of the recon- 
structed top and antitop quarks, Mf and Mf, respectively. For simplicity we neglect the 
width of the W bosons and combinatorial background in the reconstruction. The combi- 
natorial background can be estimated from Monte-Carlo simulations. We demonstrate that 
for moderate top invariant mass cuts 

\M^-mt\ < AMi (4) 

with AMt ~ 15 — 35 GeV, the phase space matching conditions are dominated by the 
NRQCD phase space contributions, i.e. they can be computed from the difference between 
the (potentially) divergent NRQCD phase space integrations without any cuts and the ones 
with the cuts in Eq. (j4]) being imposed. This is because using the MS scheme in NRQCD 
diagrams involving the unstable top propagator of Eq. ([3]) largely overestimates the contribu- 
tions from unphysical phase space regions that are parametrically away from the potential, 
soft and ultrasoft regions that can be described by NRQCD. ^ Thus the main numerical ef- 
fect of the phase space matching procedure is to remove these unphysical contributions and 
the phase space matching procedure can be carried out within NRQCD itself. The remain- 
ing hard contributions, which require the evaluation of multi-leg full theory diagrams, are 
smaller than 5 fb and can be neglected in view of the expected experimental precision. This 
simplifies the computations substantially and makes the determination of higher order QCD 
corrections feasible. Since one can assume that the situation is similar for the threshold pro- 
duction of other heavy unstable colored particles in new physics models, it is straightforward 
to generalize our results for such processes, also for the hadron collider environment. 



prescription- free cross section definitions for unstable particle production do no exist. The numerical 
impact of the selection prescriptions can, however, be frequently neglected if the width of the involved 
particles is much smaller than other relevant kinematic scales. 
^ A similar feature arises also for the computation matching conditions within NRQCD for a stable quark, 
see Ref. 



28[ for a comparison of a cutoff scheme with MS. 
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An important conceptual aspect of the invariant mass cuts defined in Eq. (jlj) is that 
aheady for moderate cuts AMt ~ 15 — 35 GeV the cut Ap on the nonrelativistic (anti)top 
three-momentum p is Ap ~ ^/2/n^[AM^ ~ 100 GeV, and thus represents a hard scale of 
the order rrit. This justifies the implementation of the phase space effects into the matching 
conditions of the Wilson coefficients. From a technical point of view, the method of compu- 
tating the phase space matching conditions within NRQCD is quite similar to using a cutoff 
scheme to regularize ultraviolet (UV) divergences in loop diagrams for the renormalization. 
The procedure therefore leads to power-counting breaking effects which could spoil the non- 
relativistic expansion. We find that this is not the case partly because the momentum cutoff 
Ap is still sufficiently smaller than the top mass rrit. In analogy to the usual renormalization 
and matching computations beyond the one-loop level, the phase space matching procedure 
also requires the determination of phase space matching coefficients of subdiagrams at higher 
orders in the loop expansion. Many technical details of the calculations carried out in this 



work are given in Ref. [29 



The outline of this paper is as follows: In Section [IT] we review the NRQCD framework 
which we use to describe the top-antitop resonance region, and collect previous results on 
electroweak and finite lifetime effects in the cross section. In Section IIII Al we discuss the 
main concepts of the NRQCD phase space matching procedure with special emphasis on the 
definition of the inclusive cross section. In Section IIIIBI we compute the imaginary type-2 
matching conditions for the Wilson coefficients of the {e~^e~){e^e~) forward scattering op- 
erators with = and a cut AMf on the reconstructed invariant masses of the top and 
antitop quarks, including also ^) relativistic corrections. The structure of the nonrela- 
tivistic expansion of the phase space matching contributions for the inclusive NRQCD cross 
section is examined in Sec. IIII CI by a comparison with the full Standard Model tree-level 
predictions for the processes e'^e~ ti ^ bbW^W~ and e~^e~ — )■ bbW~^W~ obtained from 
Madgraph jsoj. The 0{as) QCD corrections to the phase space matching conditions arising 
from ultrasoft gluon exchange and potential interactions are computed in Sec. IIVI The cal- 
culations allow to determine the complete set of NNLL phase space matching contributions 
and of a part of the N^LL corrections. The results are discussed in detail and analyzed 
specifically with respect to the convergence of the a^-expansion, and we also discuss non- 
perturbative effects. Details on the computation of the 0{as) ultrasoft corrections to the 
phase space matching coefficients are relegated to Appendix |Xl Finally, in Sec. El we analyze 
the numerical impact of the phase space matching contributions to the inclusive NRQCD 
cross section and compare their size to the other types of electroweak and finite lifetime 
corrections. Readers only interested in the main concepts of the phase space matching and 
the numerics of the final results might jump from here directly to Sec. IIII Al and then to 
Sec.lYl 



II. NRQCD FORMALISM AND PREVIOUS RESULTS 



For the work in this paper we use the vNRQCD effective field theory formalism |3lH33| to 
describe the nonrelativistic top-antitop dynamics relevant for the inclusive threshold cross 
section. In this section we briefly review the basic vNRQCD notation and ingredients with 
special emphasis on electroweak and finite lifetime effects in the cross section. We also outline 
the previous results concerning the type-1 imaginary contributions to the vNRQCD Wilson 



coefficients determined in Ref. [27|. We note that the results are general and can in principle 



also be implemented with minimal notational modifications within other formalisms such as 
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pNRQCD 0. 



Basics: 

The vNRQCD Lagrangian contains heavy quark bihnear, potential, soft and ultrasoft oper- 
ator terms. Up to the NNLL order the bihnear terms read 

+ (^p ^ Xp) , (5) 

where the fields ipp and Xp destroy top and antitop quarks with label momentum p, and Tt 
is the top quark width defined at the top quark pole. The velocity counting of the (anti)top 
quark fields for d = 4 is ^/'p ~ Xp ~ v^^"^. The term cx: is the top lifetime dilation correc- 

tion. We also included the ultrasoft gauge covariant derivative = {D^, — D) = d'^ + igA'^ 
where is the ultrasoft gauge field. The dependence on Z)^ is directly tied to the (anti)top 
three-momentum label p to all orders of perturbation theory through reparametrization in- 
variance 35|. We use the w-counting ~ rritv'^ ~ Tt, and the ultrasoft interactions start 
contributing in matrix elements at N^LL order. The term 6mt is a residual mass term which 



is of order v"^ in a short-distance threshold mass scheme The LL order terms in Eq. (|5]) 
lead to the top/antitop propagator in Eq. (|3]). 

The leading order potential term contains the well known Coulomb interaction, 



Cpot = - J2 fpl ^l'^pxlp'X~p , (6) 
p,p' ^ ' 

where 

V^'\u) = -47rCFas{mtu) (7) 

is the Coulomb Wilson coefficient for a color singlet heavy quark pair and u is the vNRQCD 
velocity renormalization scale. The velocity renormalization scale u is introduced to conve- 
niently describe the correlated renormalization group evolution of soft (~ rritv) and ultrasoft 
(~ rritv'^) effects, and its natural scaling to sum large logarithmic terms is ~ f ~ a^- The 
evolution of the Coulomb Wilson coefficient differs from the running of the strong coupling 
starting at NNLL order due to ultrasoft gluon corrections. There are also radiative correc- 
tions to the Coulomb interaction arising from soft gluon loops. In vNRQCD they contribute 
through matrix elements involving T-products of soft gluon operators. These QCD cor- 



rections to the Coulomb potential are known at 0{as) [36[ and 0{aj.) |37H39|. and have 



recently been determined even at [iMll. A discussion on the ti^-suppressed vN- 



RQCD potentials is given in Refs. [32|, |33| . We note that there are also potentials generated 
by the electroweak interactions. The dominant one is the QED contribution to the Coulomb 
potential which can be easily implemented by writing V'f\v) = — Air Cpct s{i^ti^) — 4:7r a (mtiy) . 
According to Eq. ([T]) this QED effect contributes at NLL order. There are no 0{aas) QED 
NNLL corrections that contribute to the Coulomb potential. There are also higher order 



electroweak potentials initiated e.g. by the exchange of the Higgs jo], 4o|, 41| or the Z-boson 
Eqs. (|5]) and (Q yield the LL nonrelativistic top-antitop dynamics. For predictions of the 
inclusive cross section one important ingredient is the zero-distance Green function. At LL 
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order and using dimensional regularization with c? = 4 — 2e, the unrenormahzed zero-distance 
Green function has the form (a = Cpas) 



2 





In ( 


^ —iv\ 
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hln2 + 7£; + ^ 1 

2 ^ it.^ ^\ 2^ 



m? a 1 

+ ^7T' (8) 



where 



a/s - 2(m f + 5mt) + zFf 
mt ' 

= ci/(ixlnr(x) is the digamma function, and -y/i is the e^e~ cm. energy. The divergent 
1 /e term in Eq. (|8]) is a UV-divergence in the real part of the LL order Green function and 
does not play any role for the production of stable particles. 

Top pair production currents: 

Top-antitop pair production and annihilation is described by currents involving the top and 
antitop fields. For cross section predictions up to NNLL order one needs the leading and 
subleading '^5'i currents (9p and C'p^2' respectively, and the ^Pi current Op 3, which is_p/rn+- 



suppressed compared to the leading S-wave current. These currents have the form |l7l. II8 

Oi,s = ^/A^',^■p]{^^')x%■ (10) 

To ensure electroweak gauge invariance at subleading order it is necessary to include the 
initial electron and positron fields, which leads to the ti production operators 

Ov.,p,. = [ e+ 7, e- ] O^p,. , Oa,p,. = [ e+ 7, 75 e_ ] O^,. , (11) 

where the index j = 1,2,3 is summed and the index a = 1,2,3 distinguishes between the 
different currents. The current Ov,p,i contributes at LL order for the inclusive cross section, 
and Cy,p,2 and Ov,p,3 contribute at NNLL order. Because the effective theory is constructed 
such that it describes tt production only in the threshold region and in the c. m. frame, only 
those initial e+e" states^ a^|(k')a^ f (k)|0) are allowed that fulfill s = {k + k'Y ^ Ami 
k = — k'. For simplicity we assume electron and positron to travel along the 2;-direction, 
therefore the explicit form of their 4-momenta is 



Here, ar(k) and a^,(k') are operators for the annihilation of an electron and a positron with spin r, t' 
and S-momentum k, k', respectively, and |0) is the vacuum state. 
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being the unit vector in 2;-direction. The fields e_ and e+ in Eqs. (fTT]) are now defined as 



e_(x) = ^a,(k)n,(k)e-^'=-", e+{x) = ^ <t(k') ,;^(k') g^fc'-x ^ ^^3^ 

where Mr(k) and Vr(k') denote Dirac spinors for electron and positron, respectively, and the 
momenta k and k' refer to Eqs. ( |T2l) . The sum over the c. m. energy ^/s is restricted to the 
threshold region. For simplicity we do not sum over the angles of the electron and positron 
momenta. The phase factors in Eqs. (fT3!) are defined such that they describe only the tt 
low-energy fiuctuations, 

= (t - T • = (t - ■ 

The dependence on 3-momenta in Eq. (fT4|) vanishes after the operators Ov/A.p,a have been 
applied to the initial e~^e~ state. The operators for ti annihilation are obtained from Ov/A,p,a 
by Hermitian conjugation. 

Due to the dependence of the intermediate photon and Z boson propagator on the cm. 
energy in the process e~^e~ — )■ 7*, Z* — )■ ti, we also introduce the ti production operators 

O'vli = [ e+ 7. {E/mt) e_ ] O^,, , 0%^, = [ e+ 7, 75 (E/m,) e_ ] O^, . (15) 

Here, E denotes the operator E = ido acting on the fields to the right and to the left and 
thus picks up the kinetic energy E = y/s — 2mt ~ rritv'^ from the initial e~^e~ state. So these 
operators contribute a NNLL order. Concerning QCD effects, the operators Oyj^^ ^ ^ have 
the same matching conditions and renormalization group evolution as Cy/A,p,i of Eqs. ( ITTj) 
but their Wilson coefficients differ in the electroweak contributions. Similar additional op- 
erators related to Ov/A,p,2 and Ov/A,p,3 do not need to be introduced since they would give 
contributions beyond N^LL order. 

The contribution of the currents to the Lagrangian reads 



£ = V 

•'—cur / ^ 



Cv,i Ov,p,i + Ca,i Oa,p,i + Cy^i C^^^p^i + ^a,p,i 



+ Cy,2 Oy,p,2 + Ca,2 Oa,p^2 + Cy,3 (^y,p,3 + C^,3 (^A,p,3 J + H. C. . (16) 

If QED radiative corrections are neglected, the electron and positron fields in the current 
operators act like classic fields and do not contribute to the nonrelativistic ti dynamics. 
The Hermitian conjugation (H. c.) acts on the operators in the usual way and on possible 
CP-violating phases in the Wilson coefficients as complex conjugation. It does, however, not 
act on imaginary contributions in the Wilson coefficients that are related to finite-lifetime 
contributions. Since we do not consider CP-phases in this work, the Wilson coefficients of 
operators and their conjugated counterparts are equal. At NNLL order the currents run 
only due to QCD effects. Their matching conditions, which we will generally parameterize 
at the hard matching velocity scale z/ = 1, contain electroweak contributions already at LL 
order. Including QCD, electroweak and finite lifetime effects up to NNLL order the Wilson 
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coefficients can be written in the form 



Cv/A,2{.^) 



^(l),Born / N 



(17) 



In Eqs. f[T7j) all imaginary contributions are indicated by an explicit factor of the imaginary i. 
The terms Ciiv) parameterize the QCD evolution and the hard matching conditions. Their 
explicit form at NLL order can be found in Refs. jss], H^, |43||. Note that for the NNLL 
running of Ci(i^) currently only the non-mixing contributions are fully known. Partial 

For our 



results for the NNLL mixing contributions have been determined in Refs. [20|, 1441 
numerical analysis of the phase space matching contributions to the inclusive NRQCD cross 
section we need the NNLL order corrections to the matching coefficient C\{y = 1), 



Ci[U 



1 + + + 



(2) 



(18) 



where 

,(2) 



2Cf 



TX 



In 2 



31 
24 



CaCi 



In 2 



+ 



2g? 



TT 



(19) 



and the constant k was given in Refs. j45|, |46|. The last term in Eq. ( IT9|) is the one-loop 



QED matching correction to c\(v = 1), which contributes at NNLL order according to 
the power-counting aqed ~ a^- In Eqs. (fTTl) the terms Cy°™^ denote electroweak matching 
contributions from the tree level amplitude for the e~^e~ 'j, Z ^ ti production process 
and CyJ^\ refers to the hard one-loop electroweak corrections. The results for Cy'/^i are 



23l . |24| | ) . The imaginary terms with 



elaborate and can be found in Ref. |22| (see also Refs. 
^v)a 1 contain finite lifetime effects arising from the interference of the dominant double 
resonant process e~^e~ ti ^ W~^W~bb with single resonant processes leading to the same 
final state, but having only one top or one antitop at the intermediate stage, see Fig. [H The 
terms CyJ^^ and iCyjj^ ^ are obtained from one-loop full theory diagrams and contribute at 
NNLL order in the counting scheme of Eq. ([1]) due to an additional factor of a ~ f ^. Thus at 
the order we are working they are only accounted for in the leading order current Oa/v,p,i- 
The terms iSci{A) and i5c™*(A) indicate imaginary matching contributions related to the 
experimental selection cuts generically denoted by the argument A. As we show in Sec. IIVBI 
they are required to account for a phase space matching correction of ti vertex subdiagrams 
that arise at 0{as)- For inclusive selection cuts such as the invariant mass prescription of 
Eq. (jl]) this term contributes at N^LL order. To shorten the notation we frequently drop 
the matching scale {u = 1) dependence of the electroweak matching conditions displayed 
in Eqs. f fT7|) . It is implied that the electroweak couplings are evaluated at the same hard 
matching scale as the QCD matching conditions. 
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The electroweak matching contributions from the tree level amplitude of the e'^e — )■ 
^, Z ^ tt production process read [23] (a = aqcd("^t)) 

Qt VeVt 1 ^Born . _ _ ^e^t 
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(20) 



where 



a/ 



the symbol Q/ is the electric charge, and ^3 is the third component of the weak isospin of 
the fermion /. The abbreviations s^, and denote the sine and cosine of the weak mixing 
angle, respectively. The coefficients C*y/*jf ™" arise from the 0{v'^) terms in the expansion of 
the photon and Z boson propagators near threshold using s = 4m^(l + E/mt + . . .), where 
the dots represent terms of order E"^ and higher. An alternative approach is to keep the 
exact relativistic form for the photon and Z propagators. This leads to the expressions 
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(21) 



For our numerical examinations we use these alternative definitions for C^°™. and Cyj^i^"^ 
{i = 1, 2, 3) unless noted otherwise. 

The imaginary interference coefficients Cyj^ ^ are determined from the bW cuts in the 
one-loop electroweak corrections to the e+e~ — )■ ti amplitude with on-shell stable external 
(anti)top quarks. They also contain the imaginary contribution of the (anti)top wave func- 
tion renormalization Z-factor, and this term is essential to maintain gauge-invariance. They 
have the form |21 
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(22) 
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1 + , n In 



4s?„ / \ X 



12mfsljX{4:c1 - + x) 
Qt{l - x){l + 2x){l + X + x^) 

^(1 + 12x + ^ 2x3) + -^(2 + 41x + 28x2 - x^ + 2x^) 



(23) 



where x = M^/m^. 

Forward scattering operators: 

The {e^e~){e~^e~) forward scattering operators are required to renormahze the phase space 
divergences. Their evolution accounts for the summation of Inv terms that arise specifically 
from finite-lifetime effects. The matching conditions for their Wilson coefficients account 
for selection prescriptions applied on the observed final states to define the inclusive cross 
section. These phase space matching conditions depend on the cm. energy, and we there- 
fore need to define a set of operators capable to reproduce this energy- dependence within 
NRQCD. The forward scattering operators are defined as 

of = - [ e_ 75 e+ ] [ e+ 7^ 75 {E/m,r e. ] , (24) 

and thus the O^'jj^ pick up n powers of the tt kinetic energy E = ^/s — 2mt from the initial 
e"'"e~ state. The normalization of the electron and positron fields ensures that we have 



1 J] ( 1 a.(k)<,(k') ai;;;, <t(k') 4(k) 



for the spin-averaged forward scattering amplitude. 

The contribution of the forward scattering operators to the Lagrangian reads 

Asc = + (25) 

n 

where the (^^^^(A, i/) are the Wilson coefficients. They depend on the renormalization 
velocity scale and they have a dependence on the selection cuts, generically denoted by A. 
Frequently we will use the shorter notation Cy/A = for the coefficients of the dominant 
energy-independent forward scattering operators Oy/A- For the examinations in this work 
we consider the operators Oy/A and 

QCD factorization formula: 



For the inclusive cross section of tt production close to threshold accounting for the phase 
space matching contributions up to N^LL order and for the QCD and other electroweak and 
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finite lifetime effects at NNLL order we liave tlie factorization formula 18, 22, 27 



cXinci(A) = -U^lm 
s 



n=0 

The spin-averaged lepton tensor reads 

= iE [^r'(k')7'(75)n.(k)] [n.(k)7'(75)t;.'(k') 

T,r' 

1 



(26) 



(A; + k'f (5"= 



(27) 



with the definitions of electron/positron momenta given in Eqs. (fT2l) . The quantities Al'^ 
are time-ordered products of the ti production and annihilation currents defined in Eq. (nil|l . 
Note that the electron and positron field operators, from which the operators in Eqs. ffTTl) . 
f|T5|l and flM|) are composed, only pick out the initial and final e^e~ states and do not affect 
the correlators A!"^ in any way. The explicit expressions for the A!"^ are 
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X] 



(28) 



where g = (a/s — 2mt,0). For the f ^-suppressed electroweak and finite lifetime matching 
coefficients Cy"'^^ and iCy^^ contained in Cy/^i i it is sufficient to use the LL current 
correlator All terms in the second, third and fourth lines of Eq. fl26|) are u^-suppressed 

and therefore contribute at NNLL order. This suppression originates from factors E/mt, 
p'^/mf appearing in ^i-nd two factors of cr ■ p/rrit appearing in O^^, respectively. The 
terms in the first four lines also appear in pure QCD for stable heavy quarks see e.g. Refs. [13- 
19|. The fifth line contains the phase space matching corrections related to the {e~^e~){e~^e~) 
forward scattering operators. As shown in Sec. IIII Bl the coefficients Cyj^ = Cy/A start 
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contributing at NLL, and the terms E/nitCy^^ 



are N^LL corrections. We note that all 



terms shown in Eq. fl26l) are understood as finite MS-renormalized quantities. 

We can write Af^ = 6''^/3Ai after tracing the sigma matrices of the currents in 3 dimen- 
sions. The correlators Ai can then be expressed in terms of contributions to the zero-distance 
S'-wave and P-wave Green functions of the two-body Schrodinger equation: 

Ai{v,mt,iy) = GN.^C'ia^v.mt.iy) + (vi'\iy) + 2Vi'\iy)'^ G'{a,v,mt,iy) 

[a, V, rrit, u) 

+ G'^'"(a, V, mt, v) + G'^'\a, v, rut, v) 
G^{a, V, nit, ^) ■ 



A3{v,mt,u) 



mt 



(29) 



Here, the terms [i = 2,s,r,k) are the Wilson coefficients of the f ^-suppressed po- 

tentials jssi . The correlator A2 can be related to Ai by the heavy quark equation of motion 
giving A2{v, mt, v) = v'^Ai{v, mt, f). Thus only the LL terms in Ai are necessary to obtain 
the NNLL order contributions of A2- The function G'^ is the Coulomb Green function. The 
LL approximation for G"^ is known analytically and has been displayed in Eq. ([8]). At NLL 
and NNLL order, related to the 0{as) j36[ and {stI-IsqI corrections to the Coulomb 



potential, we use the numerical results obtained in Refs. [17|, |l8|. They are based on an 
exact solution of the corresponding Schrodinger equation using computational techniques 
developed in Refs. 47, 48[. For the analytic formula for G^, see Ref. [l8|. All the relativistic 
corrections to the Green function, Q^^"^'^^^™^ are available in analytic form, see Refs. 0, 18 . 
They are computed from insertions of the f ^-suppressed potentials and the kinetic energy 
corrections. The Green function correction G"^^^ arises from an insertion of the lifetime dila- 
tion correction to the bilinear quark field operators shown in Eq. ([5]). The expressions for 
G'^'^ reads El 



G 



dil 



2mt 



^ V d d 
2 dv da 



G^{a, V, mt, v) 



(30) 



Phase space divergences and renormalization group evolution: 

Using the unrenormalized current correlators ^•'^ in the factorization formula ( l26l) leads 
to the ultraviolet 1/e phase space divergences either from insertions of f ^-suppressed op- 
erators or insertions of Wilson coefficient corrections describing finite lifetime corrections. 
The divergences are absorbed by the counterterms associated to the (e"'"e~)(e"'"e~) forward 
scattering operators Oy/A given in Eq. flM|) and are treated with the usual renormalization 
techniques known from effective theories. However, it is a novel feature that the phase space 
divergences and the anomalous dimension of the operators Oy/A are purely imaginary. In the 
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MS scheme the NNLL counterterms of the renormahzed OyjA operators have the form 27 
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Solving the resulting renormalization group equations for the Cy/A one obtains 
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where 
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2 

as{mt) 



ln^2 + ln21n(z) - Lia ('- 



(32) 



Here we have introduced the A-dependent hard scale {v = 1) matching conditions 
Cy/yi(A, 1).^ They are determined by the phase space matching procedure as described in 
the following sections and incorporate the information on the experimental selection cuts as 
well as the contributions from background diagrams. The phase space logarithms resummed 
in Eq. f l3ip correspond to logarithmic terms involving ratios of the hard scales rrit, A and 
nonrelativistic kinematic scales. They contribute at order ~ in the inclusive cross 
section. So compared to the LL cross section, which counts as a'^v ~ v^, the phase space 
logs constitute NLL contributions. This is expected since the phase space divergences arise 
from matrix elements contributing at NNLL order. In the following sections we determine 
the matching conditions Cv/a{A., 1), which contain the details of the selection cuts that are 
applied for the definition of the cross section. As explained in the introduction, for inclusive 
cross section definitions the information on these selection cuts represents hard effects within 
the NRQCD framework. 



III. CONCEPTS OF PHASE SPACE MATCHING 

In this section we discuss the main concepts that go into the phase space matching on 
the basis of an explicit computation of the phase space matching conditions for a cut on the 



^ Frequently wc drop the argument V = 1' in the coefficients Cy/^(A, f ) to simplify the notation. 
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invariant masses of the reconstructed top and antitop quarks. For simplicity we neglect in 
this section the effects of the strong interactions, i.e. we set = 0. A number of additional 
issues arise if QCD effects are included. QCD corrections to the phase space matching 
conditions are computed and discussed in Sec. HVl 



A. Basic Setup 

We define what we call the invariant masses of the top and antitop quarks through the 
reconstructed masses of the bW~^ and bW~ systems coming from the top and antitop decays: 

= p1 = (Pb + Pw+f , = pI = in + pw-f . (33) 

Without QCD effects related to jet and soft particle emission this definition is unambiguous, 
and we neglect combinatorial background in the following. The latter would have to be 
determined from Monte-Carlo studies in addition to the examinations carried out in this 
work. We consider selection cuts on the top and antitop invariant masses of the form 

(mt - AMJ < Mtj < {nit + AMt) . (34) 

We note that to keep Eq. ambiguity-free once QCD effects are accounted for, it is 
necessary to employ a short- distance top quark mass definition for rrit that is suitable for 
reconstruction. Such mass definitions have a low infrared cutoff scale i? ~ to avoid the 



pole mass (9(Aqcd) renormalon [15|. An example is the jet mass introduced in Refs. [49|, |50 
The MS mass does not belong to this class of short- distance masses since it has i? ~ m^. 

The constraint on the top and antitop invariant masses f can be translated into a 
condition on the off-shellness p'^^ — mf appearing in the top and antitop propagators, see 
Eq. ([3]). At LL order in the nonrelativistic expansion the constraint on the nonrelativistic 
off-shellness in the NRQCD propagator has the form 

Pit 

2mt 

The relativistic NNLL order corrections to these constraints are given in the second part of 
Sec. IIIIBi Through momentum conservation Eq. (1351) leads to a constraint on the phase 
space integrations for the cross section as illustrated in Fig. |2l From Eq. (l35l) we see that 
the ultrasoft and the soft phase space momentum integrations are limited by the scales AMt 
and \/2mtAMt, respectively. In this work we consider moderate invariant mass cuts with 
AMt ~ 15 - 35 GeV. Thus we have 



,2 



mtir < AMt -15-35 GeV , rritv < ^2mtAMt ~ 70 - 110 GeV , (36) 

and the upper bounds of integration are substantially above the generic w-scaling of the soft 
and ultrasoft momentum components. Since ^j2mtAMt is parametrically of order irtt we 
use for our bookkeeping the counting y/2mtAMt ~ m^. This also implies that AMt ~ ^t- 
Within the vNRQCD framework this counting scheme is natural since, due to the pull-up 



mechanism [33|, |51[, the ultrasoft scale is directly connected to the hard matching scale via 
RG evolution without an additional matching at the soft scale. In this counting scheme the 
phase space constraints are incorporated through the NRQCD Wilson coefficients. On the 
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FIG. 2: Phase space integration region spanned in the = (poi p) variables for tt pair production. 
In terms of and the external energy E = ^/s — 2mt the top and antitop momenta read p^^ = 
{nrit + E/2,0) ziz pf^. The red solid lines correspond to the top and antitop nonrelativistic on- 
shell conditions; in the limit = the phase-space shrinks to the point where the two on-shell 
lines intersect, {po, |p|) = (0, ^/mE). Gray, light-gray and white areas correspond to the double- 
resonant, single-resonant and non-resonant regions, respectively. The gray area represents the 
phase space region compatible with the invariant mass constraints. We have chosen E = 5 GeV 
and AMt = 20 GeV for this picture. 




FIG. 3: Graphical illustration of the factorization formula (I26p for as = 0. The tt production 
operators Oy^^^ ^ ^ for a = 1,2,3 are represented by a crossed circle, gray circle and gray box, re- 
spectively. Diagrams in the second line correspond to the contributions from the forward scattering 
operators. 



other hand, numerically the scales AMt and ^/2Jr^AM^ are sufficiently below the top mass 
scale such that all tt phase space configurations that pass the invariant mass constraint can 
still be adequately described by NRQCD. This fact is crucial for the phase space matching 
method we describe in the following. 

When QCD effects are neglected the factorization formula in Eq. fl26|) can be illustrated 
graphically as in Fig. Ei For = and in the absence of potentials the phase space matching 
conditions can only contribute to the Wilson coefficients Cyj'^ of the {e~^e~){e~^e~) forward 

scattering operators Oy/j^- This is because off-shell (anti)top phase space contributions 
(corresponding to the light-gray and white areas in Fig. [2]) do not match to the operator 
structure of the {e~^e~)(tt) currents. From Eq. ( !26|l and Fig. [2] we see that for the determina- 
tion of the Wilson coefficients Ci^L we need to know the result for the inclusive cross section 
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(a) (6) (c) (d) 

FIG. 4: (a) Full theory diagram for e^e^ — t- — t- ^bl^"^!^^. (b,c) Typical single-resonant full 
theory diagram for e~^e~ — t- tbW~ — ?• bhW^W" . (d) Typical full theory diagram for e^e~ — )• 
bbW~^W~ without top or antitop quarks as intermediate states. 



'^hid^i^) with invariant mass constraints. In the common approach to matching compu- 
tations CT^^^^i^) is computed in the full relativistic theory. After the result is expanded 
nonrelativistically using the counting of Eq. ([T]) one can identify the pieces belonging to the 
Wilson coefficients Cy^^. On the other hand, as mentioned above, the ti phase space regions 
passing the invariant mass cuts can be determined within the nonrelativistic expansion. We 
therefore write the expression for the inclusive cross section as a sum of two terms, 

Ccr°(A) = a-=Q°cD(A) + <enr°(A) . (37) 

Here, ctnrqcd cross section computed from NRQCD Feynman rules with the (anti)top 

invariant mass constraints being applied for the phase space integration. The parameter A 
is related to the invariant mass cut AMt and we use the formal counting A rrit according 
to the discussion above. The exact definition of A will be discussed below. In this compu- 
tation the {e~^e~){e~^e~) forward scattering operators do not contribute, and the resulting 
expressions are just the nonrelativistic expansions of full theory squared matrix elements 
containing the square of the double resonant diagram e~^e~ ti bbW~^W~ (see Fig. H^) 
and the interference of the double resonant diagram with the diagrams for e~^e~ — )■ bbW^W~ 
having only either the top or the antitop in intermediate stages (see Fig. IDd and c for typ- 
ical diagrams). As we show in Sec. IIII Bt the contributions to the Wilson coefficients Cy^^ 
that result from '7nrqcd(A) are local (i.e. energy-independent) and only depend on powers 
of Tt/rrit and A/rrit. While Tt/rrit ~ w^, which obeys the natural NRQCD counting, the 
Ajvut term is of order unity and can - as we show in Sec. IIII Bl - lead to power counting 
breaking contributions for insertions of operators that are higher order in the nonrelativis- 
tic expansion. However, we find that the numerical effects of the power-counting breaking 
contributions are very small and do not spoil the nonrelativistic expansion. This is partly 
due to the fact that the phase space cutoff A is sufficiently smaller than the convergence 
radius of the nonrelativistic expansion. We refer to this feature as "mild" power-counting 
breaking. 

The remainder contribution of the inclusive cross section, (T^=jj7°(A) accounts for all other 
contributions to the full theory matrix element. This includes for example pure background 
t^e~ — bbW~^W~ diagrams, see Fig. HJi for a typical diagram, and also the square of the 
single-top diagrams in Figs. and c. In Sec. IIII Cl we determine the remainder contribution 



from a numerical analysis using MadEvent |30|. We demonstrate that the remainder con 



tribution is very small and can be neglected in view of the experimental precision expected 
at a future linear collider (see Sec. H]). Restricting the MadEvent amplitude to the diagrams 
with an intermediate tt pair we also show the excellent approximation that is provided by 
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the nonrelativistic computations in o"^rq(^j-,. 
B. NRQCD Phase Space Matching 

In this section we compute o"nrqcd' NRQCD cross section for a,; = with a cut AM^ 
on the invariant masses M^ j, see Eq. flMl) . As explained above, we treat AM^ and A as a 
hard scale. 

Leading order diagram: 

We start with the leading order NRQCD diagram to set up the notation and explain our 
method of computation. Technically, the least involved method to determine the form of the 
phase space integral and the proper normalization factors is to use the factorization theorem 
of Eq. ( I26l) and apply the cutting rules on the (anti)top propagators in the NRQCD current- 
current correlators. After identifying the top and antitop momenta according to Eq. (155]) it 
is then straightforward to derive the expression for the phase space integral and the phase 
space boundaries compatible with the invariant mass constraints. 
We start from the form of the leading order current correlator 

>l?'"^=°(w,mt,z/) = 6N,G\a = 0,v,mt,u) 



It is easy to identify p^^ = {rrit + E/2,0) ± p^^ as the top and antitop four-momenta. Using 
the cutting rule for the unstable (anti)top propagators 



2 ±^0-^+^^ 



-2Im 



(39) 



and recalling the form of the invariant mass constraints in the nonrelativistic limit given in 
Eq. (135|) we obtain 

d^P 



a°'srA(A) =Ar,((c^r)' + (ctr)') 



A(A) ' y2^P0 2mt ^ 2 ) \2 2mt ^2 



^2" ^0 2^ + ^Tj(^2" 2^ ^Y 



where A (A) stands for the phase space constraint 

E p2 



2 ^ 2mt 



< AMt (41) 



as derived before in Eq. ( I35l) . The form of the allowed region in the (po;P)-plane is shown 
in Fig. |2] and a graphical illustration of the computation in Eq. f HOj) is depicted in Fig. |5l 
The expression in Eq. P0|) agrees with the result of a full theory cross section computation 
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'A(A) 



FIG. 5: Leading order diagram for the phase space matching computation. The red dashed hne 
in the diagram on the LHS means that we cut through the (unstable) top Hues and integrate 
over the phase space region defined by A (A). This is equivalent to the full theory cross section 
computation from the tree-level matrix element for e^e~ tt ^ bbW~^W~ (shown on the RHS) 
in the nonrelativistic limit. The average over the electron and positron spin states is implicit. 

from the tree-level matrix element for e~^e~ tt ^ hhW^W~ in the nonrelativistic limit. 
In this computation one has to first carry out the hW^ and hW~ phase space integrals as 
a function of the top and antitop invariant masses, respectively. The proper nonrelativistic 
limit is obtained by setting the invariant masses of the hW^ and hW~ systems to (which 
gives two factors of the on-shell width Fj) and by taking the nonrelativistic limit of all other 
remaining terms. Upon contraction of all indices in Dirac space and resumming the width 
terms into the top propagators - which is required from the counting in the double resonant 
kinematic region - one arrives at Eq. (140|) . 

A more compact representation of the ti phase space integral in Eq. ( 140|) is obtained by 
switching variables to the nonrelativistic invariant mass variables ti and t2 defined by 

ti,2 = 2mJ^±Po-T^). (42) 



Inverting the relations gives 



2 ' 2mt 



tl — t2 2 7-1 ^1+^2 



= P'-Em,~- (43) 



Using the Jacobian dp^d^p = n/ {2mt) x y mtE — |(ti + ^2) dtidt2 the expression in Eq. ( 140 
takes the form 



^3p2 r \ nitE - Uti + 12) 



A(A) 

(44) 



with 



A%t„t2) = 1 for a°'-=°D. (45) 
The integration region in (ti, t2)-space has the form 

A(A) = |(ti,t2) eM^ : <A') A (0<mtE-^{h + t2)^Y (46) 
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FIG. 6: Allowed phase-space region in the ti — t2 plane. The left and right panels show the cases 
E < and E > 0, respectively. The triangular-shaped blue shaded regions corresponds to the exact 
relativistic phase-space region allowed by the kinematics and the invariant mass cuts in Eq. (j34p . 
The orange region located on the upper right of each panel is forbidden since |p| must be real. 
The area inside the dashed lines corresponds to the points that pass the conditions (|46|) . When 
the relativistic corrections to the cut A are included according to Eq. (j52|) . the allowed region is 
the one indicated by the dotted lines. 



where = 2mtAMt. In Fig. [6] the allowed region in the ti-t2 plane is illustrated by the 
area within the dotted lines for E < (left panel) and E > (right panel). The second 
condition in Eq. f HB]) for ti 2 is from kinematics to ensure real values of the momentum 
IpI = rritE — {ti + ^2)72. Beyond the nonrelativistic approximation the relation between 
A^ and AMf receives additional relativistic corrections which are discussed below. We have 
to evaluate Eq. (jUl) in a nonrelativistic expansion for 

iTitE , nitTt ~ m^v"^ < A^ . (47) 

To this end, we carry out an asymptotic expansion based on the four regions (^1,^2) ~ 
{mfv'^,m^v'^) (double resonant), (^1,^2) ~ {mfv'^,A'^) and (^1,^2) ~ {A'^,mfv'^) (single reso- 
nant), and (^1,^2) ~ (A^,A^) (hard). We obtain 

?/ , 2^2 Fi 4 + 2y2sinh^Vl) mjF? 



4'rqcd(A) = 2Ar. ((C^-)^ + (C^D^) ^ (im M - ^ ^ + 



3 7r2 A3 



2V2 nitETt 



Ov^iri I (48) 



Stt A3 V A5 

where v = ^^{E + iTt)/mt. The first term in the parenthesis is the well known NRQCD 
Born cross section obtained from the unrestricted phase space integration. It constitutes 
the leading order cross section and is 0{v), see Eq. ([H]) for = 0. The second term in 
Eq. ( HHl) proportional to Ft/A is the dominant phase space correction and of order f^, i.e. 
it contributes at NLL order. We also see terms proportional to mtT'f/A^ and mtETt/A? 
which contribute at N3LL order. These NLL and N3LL corrections are subtracting the 
phase space contributions that do not pass the invariant mass constraints, represented by 
the white regions in Figs. [61 One may wonder why these corrections are polynomial in E 
and Tt and not a non-trivial function oi E/Tt given that this region also has single resonant 
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FIG. 7: One-loop renormalization condition for Cyj'j^{A). 



regions where either ti or t2 are resonant, i.e. of order m^v"^. This can be understood from the 
form of the integrand in Eq. ( I44l) which shows that it is impossible to generate a nontrivial 
E /Vt dependence if ti ^ t2 ~ m^t>^ or ^2 ^ ~ mfv"^. Thus the contributions that arise 
from the invariant mass constraints, and in fact all selection criteria that do not cut into the 
double resonant phase space region (^1,^2) ~ (m^f^,m^w^) represent hard contributions that 
can be matched onto the local {e^e~){e~^e~) forward scattering operators Oy^j^. Accounting 
for the form of Eq. (126|) it is then straightforward to identify the contributions to the Wilson 
coefficients Cv/a{-^, 1) and Cyjy^{A, 1). A graphical illustration of the procedure is shown in 
Fig. [71 We obtain the results 
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(49) 



and 



/=.(1),0 



(A,l) 
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1^ 



47r Svr A^ 



(50) 



The coefficients 6*°^^ (A, 1) contribute at NLL order and the Cy'j'^{K, 1) at N^LL order. The 
superscripts "0" are a reminder that QCD corrections are neglected. 



Relativistic corrections: 

We now determine the phase space matching contributions from 0{v'^) relativistic corrections 
to the NRQCD cross section. Since they come with one additional factor of p^/m^, they are 
suppressed by according to the counting p ~ nitv. However, since the imposed invariant 
mass constraints act like a hard momentum-cutoff for the NRQCD phase space integration, 
power-like dependences involving the cut parameter A appear as a consequence of employing 
the nonrelativistic expansion. This occurs whenever the nonrelativistic expansion leads to 
integrals that are divergent for |p| — j- 00. In effective field theory computations where a hard 
cutoff regularization is used, this means that once power-counting breaking occurs from the 
insertion of a higher dimensional operator O, Wilson coefficients contributing to a lower 
order than the operator O might have to be modified during the matching procedure to 
the full theory. This happens e.g. in lattice computations and, while it is admittedly not 
esthetically elegant, it only represents a technical subtlety from the field theoretic point 
of view. For the method we use for the phase space matching, however, the occurrence 
of power-counting breaking terms that are numerically large could be disastrous since our 
method relies on the assumption that the ti phase space passing the selection cuts can be 
computed reliably in the nonrelativistic expansion. This restricts the size of the invariant 
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mass cut AMt we can handle within our method for the phase space matching procedure. 
We show in Sec. IIII CI that for the invariant mass cuts AM^ between 15 and 35 GeV which 
we consider in this work power-counting breaking effects are very small and do not cause 
any problems. 

There is an important relativistic effect related to the invariant mass constraint given 
in Eq. fl4ip . These corrections modify the bounds on the ti^2 integrations in Eq. (H^ . To 
determine these corrections one needs to solve the equations 

[nitT^M^^ = 0, 

(mj±AMt)'=0, (51) 

where E = a/s — 2mt. The condition mtE — ^{ti +^2) > remains unchanged. The allowed 
phase space region in the ti — t2 plane accounting for all relativistic effects is illustrated by 
the blue shaded regions in Fig. |6l Since an exact analytic integration within the relativistic 
phase space boundary is not required for the precision aimed for in this work, we derive the 
following approximation. The main effect of accounting for the exact boundaries is a change 
in the lower boundaries of the ti^2 integrations. The approximation can be made because the 
single resonant regions where either ^2 or ti are zero {i.e. the regions close to the negative 
X- or y-axes inside the blue shaded areas in Fig. E]) give the largest contributions. For the 
determination of the new lower boundary for the ti 2 integration one can take, for example, 
the first of the relations in Eqs. f lSTjl . Switching to the t- variables and setting for example 
^2 = 0, it is straightforward to find the relation 

|ti,2| < A' = 2mtAMt - ^{AMtf - ^EAMt + ... (52) 

The associated region in the ti-t2 plane is indicated by the dotted lines in Fig. |6l The 
omitted terms in Eq. fl52|) indicate: (a) terms involving the energy E, which contribute 
beyond N'^LL order to the cross section and thus can be safely dropped for the range of 
energies E considered in our analysis, and (b) terms which give power-counting breaking 
contributions similar to those arising from the N^LL order kinetic energy corrections, that 
are shown below to be numerically negligible. Apart from the neglected terms in Eq. flS^ . 
the difference between the exact relativistic phase space region (with curved boundaries) 
and the approximation given by the dotted lines in Fig. |6l can also be shown to yield terms 
beyond N^LL, or terms of type (b). 

The other contributions from 0{v'^) relativistic corrections to o"nrqcd(^) be cast into 
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the form of Eq. (jHj), where the functions A*(ti,t2) read 



Akin,0 _ P [ i_ 



, dil,0 



, P-wave,0 



m? (t? + m?r?)(t^ + m?r?)' 

o/^Born/^Born i O/^Borri/^Born 2 

2 (c^,r)2 + {cfyy p2 



3 (C^^r)^ + (CaT)^ ^ 



2 ' 



/^Born/^int i /^Born/^int , , , 
into _ ^V,l '-^y,! ^^,1 "-"A,! 11+12 



(53) 



Here A"^™'*^ comes from the insertion of the kinetic energy correction p^/(8mj) and A*^'''" 
from the time dilation correction —i Ft p^/ (4m^) contained in the quark bihnear Lagrangian 



Eq. The term A*" ' comes from the insertions of the f -suppressed 5- wave current Op^2 
and A^-"^™'° from the P-wave current Op,3, see Eqs. ( fTOl) . Finally, the function A™*'° arises 
from interference contributions of the double resonant amplitudes e~^e~ ti ^ hhW^W~ 
with those where only either the top or the antitop appear at intermediate stages, see 
Eqs. (122|) and (123|) . Using the methods of the previous subsection to compute the corrections 
to the cross section with invariant mass cuts we obtain the following expressions for the 
contributions to the Wilson coefficients (^^/^(A, 1) and (7^j^(A, 1): 

r^nfi(K 1^ _ 9,- Ar rr-Born ^^ f 9 T^A -86 + 105^2 sinh'^(l) x 
Cy,^ (A, 1) - 2 z AT, (C^/^, J ^ + — ^ j , 

r<dii,o.A 1^ _ 9. Ar rr^Born n2 Sm^ / 2v^ T^A 4(2-3v/2sinh-\l)) x 
C^/^(A, 1) - 2^7V, (C^/^,,j T6^ I - ^ + ^^^^ ' 

r-^O.. n - 9,- /V 9rBorn ^Born ^1 T^A -2 + 3v/2 siuh"^ (1) X 

C^/^(A,1) - i^7V,2C^/^,iC^/^,2^(^— -5 — ^j, 

^p-wavc,o.. n - ^VAr rr^Born^2 ^ v/2 T.A , -2 + 3 v/2 siuh'^ (1 ) X 
1) - 3 ^ (^v^M,3j (, V + ' 

r<mt,o. . . N p , ;v 9 + -2 + 3y2sinh-^(l) m^EA 
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and 



c«*»(A.i) = 2iiv4c!?;j;0^|f^^ 



/=,(1),d2,0/ . _ r, • AT r) ^Born ^Born f^t "^V^ 

/=,(l),P-wavc,0/ . _ 4 /^Born \2 3a/2 

C'i^/r'^A, 1) = -2 z AT, 2 C^/TiCl;;^^,,^ ^ , (55) 

where Cy°™2 = "V^C'y/™!- As anticipated for the results contributing to Cy/A(A, 1), 
apart from the N'^LL terms oc T'^/mf, there are in Eq. fl5^ terms oc TtA/rrit that are 
power-counting breaking since they contribute at NLL order. These contributions featuje a 
relative factor A^/m^ = 2AMt/mt = 0.15 — 0.4 with respect to the NLL order terms in Cy/^ 
of Eq. (149|) that come from the cross section in the nonrelativistic limit. For the interference 
contributions this statement also applies because Cyjj^^ ~ Cyy^^r^/m^. The small size of 
this factor is one reason why these power-counting contributions do not spoil the quality of 
the nonrelativistic expansion. 

For visualization we show in Fig. [HI as a function of AMf, the numerical contributions 
to the inclusive cross section from the NLL phase space matching coefficients Cy^^{A, 1) + 

E /rritCyj^ {A,l) in Eqs. ( H9!) and Eqs. (150|) (red sohd line), and those from the different 

NNLL relativistic corrections, ^^^^(A, 1) + E/mtCy^'l''^{A,l) {i = kin, dil, w^, P- wave, int) 
given in Eqs. (!54|) and (155!) (blue, green, brown, magenta and cyan dashed lines, respectively). 
The difference of dashed lines between left and right panels illustrates the size of N^LL effects 
in the E /rritCyj^'^ {A, 1) terms. Note that the various contributions have different signs and 
that in Fig. [8] only their absolute value is displayed. For comparison, the solid black line 
shows the NRQCD cross section without phase space cuts. The sum of all the contributions 
in Eqs. ( 15^ and ( l55l) is displayed as the dotted black line, which almost overlaps with 
the interference contributions. We see that except for the interference contributions all 
phase space matching corrections in Eqs. (l54l) are an order of magnitude smaller than those 
obtained from the leading order NRQCD cross section in Eq. (149|) . This shows that the 
power-counting breaking terms oc TtA/mf are small and do not spoil the nonrelativistic 
expansion. Interestingly, in the sum these phase space matching corrections also cancel to 
a large extent due to their different signs. In Fig. [8] we also see that the contributions from 
the interference coefficients, Cy^'^{A, 1), are the by far largest terms that come from NNLL 
order relativistic insertions. For AMt > 35 GeV they are comparable to the contributions 
of ^^/^(A, 1). Since we have numerically that 2C|^J°Ci?* + 2C|°™C^"\ = -4.7 ((C^°i™)2 + 
(C^°™)^) Ff/mj, we see, however, that the large size of the interference terms comes from 
the size of the coefficients CyJ^ ^ and is not related to power-counting breaking effects. As 
far as the size of higher order relativistic corrections are concerned we thus also expect a 
good perturbative behavior for the interference effects. 
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FIG. 8: Absolute values of the contributions to the inclusive cross section from the leading order 
diagram phase space matching coefficients, Cy^^{A,l) + E/mtCyj^{A,l), found in Eqs. (09]), 

and (jSOp (red solid line), and those from the different NNLL relativistic corrections, (7^^^ (A, 1) + 

E /nitCy^^''^ {A,l) {i = kin, dil, v^, P- wave, int) given in Eqs. ([5^ and ([55]) (blue, green, brown, 
magenta and cyan dashed lines, respectively). For comparison we also display the corresponding 
contributions from the tree-level NRQCD cross section without phase space matching contributions, 
f^NRQCD(^ = oo) = 2iVc ((C^°i™)2 + (C^°™)2)m2/(47r)Im(iw) (black horizontal lines). The left 
panel corresponds to E = —5 while for the right we have E = 5 GeV. The values chosen for the 
input parameters can be found in Eq. (|63p . We have used the energy- independent form of the 
coefficients C^/°^° and C^y/f"''' {i = 1,2,3) given by Eqs. ([20]) . 



To examine in an example the behavior of power-counting breaking effects from operator 
insertions beyond NNLL order we now consider the phase space matching contributions aris- 
ing from the N^LL order kinetic energy corrections. These emerge either from two insertions 
of the 0{v^) kinetic energy operator p^/ {8mt)^ or from one insertion of the 0{v^) subleading 
kinetic energy operator — p^/ {IQrrit)^. The respective expressions for the functions A*(ti, ^2) 
read 



^km,0,2xp''/(8m3) 



32mj 



3tl 



{ti + m^r^y 



3tl - 5m?r2 



(ti 



^km,0,lxp6/{16m5) 



2ml 



kin.O 



(56) 



and the corresponding contributions to the Wilson coefficients Cv/a{^i 1) and C|y^(A, 1) 
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FIG. 9: Absolute values of the phase space matching contributions to the inclusive cross section 
from the NNLL (from Eqs. (j54ll55p . blue lines), and N^LL insertions of the kinetic energy operators 
(Eq. ([571 [58j) . green lines), as a function of AMt for E = —5 (left panel) and 5 GeV (right 
panel). The solid lines show the contributions from Cy/yi(A, 1) and the dashed lines show the 

contributions for C'y/^(A, 1) + E /mtCyj^{A, 1). For comparison we also display the corresponding 
contributions from the tree-level NRQCD cross section without phase space matching contributions 
(black horizontal lines), and the dominant phase space matching corrections from Eq. (149^ and 
Eq. (f50l) (red lines). We have used the energy-independent form of the coefficients Cy"^- and 

Cv/A,r = 1' 2, 3) given by Eqs. ([20]). 
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and 



53 FjA^ 



5 (754 - 483^2 smh~^(l)) F^A 



^kin,[lxpV(16mf)],0, . ^^ _ ^ . , . /^Bom\2-^/ 



11767r2 mf 

7m? / 11 F,A3 
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647r V42V27r 

514- 315y2sinh-^(l)F2A . 
+ 140^^^ ^^^^ 



:,(l),[2xpV(8m3)],0.. . _ . BornN2 63m^^l5^FiA 



In Eqs. ( 1571) and ( l58l) we have displayed all terms that contribute at NLL and N'^LL order. 
They all are proportional to powers of A/mt and power-counting breaking. In Fig. |9] the 
numerical contributions of the kinetic energy corrections to the inclusive cross section from 
the NNLL (from Eqs. (I54l [55]) . blue lines) and N'^LL insertions (Eq. (ISTf [581). green lines) of 
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the kinetic energy operators are displayed over AM^ for E = —5 and 5 GeV. The respective 
sohd hnes show the contributions to Cv/a{-^, 1) and the dashed hnes show the contributions 
to Cv/a{-^, + E /mtCy^j^{A, 1). For comparison we also display the corresponding contribu- 
tions from the tree-level NRQCD cross section without phase space matching contribution, 
^nrqcd(^ = oo) = 2iVe ((C^°r)2 + (C5°r)^) m2/(47r) Re(u) (black lines), and the dominant 
phase space matching corrections C°^^(A, 1) from Eq. f HOj) (red lines). The red dashed hnes 
account for the energy dependent terms in Eq. (150|) . that also originate from the leading 
order diagram. The results show again that despite the existence of contributions that are 
power-counting breaking the numerical impact of these terms is very small. The phase space 
matching corrections from the NNLL kinetic energy insertions never exceed the level of a 
few fb and those from the N^LL insertions, which are purely power-counting breaking at the 
order we consider, are of order 0.01 fb or smaller. Morever we see that the contributions to 
Cv/Ai^, 1) and to E/mtC^y^^{A, 1) coming from NNLL kinetic energy insertions are similar 
in size, as would be expected in the absence of power-counting breaking terms. This confirms 
that power-counting breaking effects are not of any concern. In particular we see that the 
quality of the nonrelativistic expansion is excellent and the corrections from insertions of 
higher dimensional relativistic operators are numerically compatible with the nonrelativistic 
f-counting. 



C. Full Theory Analysis 

In this section we analyze results for the inclusive cross section obtained from relativistic 
amplitudes in the full theory for = 0. In the first part we compare the NRQCD prediction 
using the factorization formula of Eq. (126|) at NNLL order with the corresponding inclusive 
cross section from the full theory accounting only for the amplitude e^e~ tt ^ bbW~^W~. 
In this case the remainder contribution a^^"^^{A) vanishes since this amplitude is fully ac- 
counted for in the usual NRQCD matrix element and matching computations. We show that 
the NRQCD factorization formula provides an excellent approximation to the full theory pre- 
diction for invariant mass cuts AMt = 15-35 GeV. We also demonstrate that the NRQCD 
prediction in the MS scheme without the phase space matching contributions^ overestimates 
by far the full theory predictions. This shows the importance of the phase space matching 
procedure and also explains why the NRQCD contributions turn out to be the dominant 
terms in the phase space matching corrections. In the second part of this section we compare 
the full Standard Model prediction for e'^e~ — )■ bbW~^W~ with invariant mass cuts with the 
NNLL NRQCD prediction including the NRQCD phase space matching contribution. We 
demonstrate that, as indicated before, the remainder cross section a^^^^{A) is small and 
only amounts to at most several femptobarn. For the precision expected for threshold cross 
section measurements at a future linear collider the remainder contributions can therefore 
be safely neglected. For the full theory computations carried out in this section we use 



MadGraph and MadEvent [30|]. At this point we note that the amplitudes generated by 
MadGraph are at the tree-level and use the fixed-width scheme for the top quark propaga- 
tor, + rrit) I {q^ — mf + irritTt). This expression for the top propagator is the correct form 



^ All NRQCD predictions for top pair threshold production that can be found in the previous literature 
were carried out in this approximation. 
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in the resonance limit and compatible with the nonrelativistic expansion of NRQCD at the 
order we are working. We note that a discussion on invariant mass cuts for the threshold 
production of a W^W~ pair in e^e~ collisions based on full theory calculations similar to 



the ones given here was presented in Ref. [52[. We also refer to Refs. |53l. |54 
Analysis for e+e~ — tt — )■ hhW^W^: 

Since the amplitude for e^e~ Z ti ^ bbW~^W~ contains a ti intermediate state, 

its contribution in the tt threshold region is fully described in NRQCD. Thus as long as 
invariant mass cuts are applied in the region where the nonrelativistic expansion is valid, 
the phase space matching contributions can be computed entirely within NRQCD and the 
remainder contribution to the cross section, (T^-'jj^°(A) is zero. Thus the comparison of the 
NRQCD prediction of the inclusive cross section with invariant mass cuts based on Eq. (126|) 
with a fully relativistic tree-level computation based on the same amplitude serves as an 
important numerical check of the nonrelativistic expansion and of the NRQCD phase space 
matching method. In this section we carry out this check for the tree-level full theory cross 
section (i.e. without accounting for one-loop electroweak corrections) and at NNLL order 
in the nonrelativistic f-expansion. We note that the amplitude for e~^e~ — )■ bbW~^W~ is 
not gauge- invariant if one includes only the diagrams with a ti intermediate state. ^ To be 
definite we therefore pick the unitary gauge for all the calculations in this section. 

For the relativistic calculations we employ MadGraph for the amplitude generation and 
MadEvent with 10^ events to numerically compute the cross section with cuts on the top and 
antitop invariant masses. For the NRQCD calculation one employs the factorization formula 
of Eq. (126!) neglecting all QCD effects from matching coefficients and potentials and the real 
one-loop electroweak corrections to the Wilson coefficients CyJ^^ of the {e~^e~){tt) current 
operators. One also needs to neglect the imaginary interference contributions iCyJj^ ^ to the 
Wilson coefficients of the {e~^e~)(tt) current operators (see Eqs. fl22|) and fl23|) ) because we 
do not account for diagrams with either only a top or an antitop at intermediate stages. 
However, there are off-shell corrections to the (anti)top decay from the full theory relativistic 
amplitude which in the matching procedure for iCyjj^ ^ are related to the imaginary part of 
the top wave-function renormalization constant in unitary gauge. These imaginary matching 



contributions have the form 27 



(zC-)^-- = -^ ^2^,4 (4^2 _ ^) [^QeQtsli^ -x + Ax^- 9x^ + 2x^) - 4Q..^(1 - x^ 
+ 3Qtsl{l + x^ - 2x'^) - 1 + x"^ 

(iC- )^-- = ^ ""'/jjfl [3g,4(l + 2x') - 1 + , (59) 
32xmfsl,[4:C^ - x) 1 J 

and need to be included. The resulting expression for the NRQCD cross section at NNLL 



2 



In the nonrelativistic expansion the gauge-dependence of the resuhs based only on the amplitude for 
e+e" ^ 7, Z ^ tt ^ bbW+W- starts at NNLL order. 
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order is remarkably simple and has the form 



NRQCD.tt - -c^vv-y,i ; - v-a,i ; Jy27T ' ^ tt^A 48^2 



-Re{iv) H -lm{iv^^ 



Stt ' ' 487r 



+ (c^,r(^^a)^""°' + cir(c^'i)^""") (^^^1?^ + ^^M^))) , (60) 

where the subscript ti is a reminder that for the amplitude only diagrams with a ti pair 
in the intermediate state are considered and where we use the s-dependent convention for 
the current coefficients Cy°™^ given in Eq. fl2T]) . The corresponding LL cross section is just 
the well known nonrelativistic lowest order expression which does not contain phase space 
matching corrections 

<mclu = {(C^rf + iC^n') ^ Im M , (61) 
and the NLL cross section reads 



2 



-nrqS(a) = iVc {{Krr + (^?r)') (l7i-(-) - 

v^r^A 



TT^A 48y27r2 



37r2 



+ ((Cyr)' + (CAr)')- 

I //^Borri/'/^int \Zt,uni , /^Born //^int \ ,uni\ ^^^^t ■'^ A ^«oN 

+ ('-y,i ('-y,i) +^A,i (W,i) ) — — )■ (°2) 

In the upper panels of Fig. [10] the full relativistic cross section a^f~^{A) obtained from 
MadEvent (red lines) and the LL (lower green dotted line), NLL (blue dash-dotted lines) 
and NNLL order (blue dashed lines) NRQCD cross sections of Eqs. fl6T]) . fl62|) and fl60|) . are 
shown for AMt = 15 (left panels) and 35 GeV (right panels). We use the following set of 
parameters for the computations: 

mt = 172.00 GeV, Tt = 1.4614 GeV, 
Mw = 80.419 GeV, Mz = 91.188 GeV, 

= Mw/Mz, a = 1/132.51. (63) 

We see that the quality of the nonrelativistic expansion is excellent. Moreover, at the scale 
of the figures no differences between the full relativistic and the NNLL order NRQCD cross 
section are visible. The corresponding lower panels show the difference of the relativistic 
and the NNLL NRQCD cross section, Aau{A) = a^f=° (A) - a^'g^^^^^^ (A) , which illustrates 
the small size of corrections coming from beyond NNLL order in the NRQCD computation. 
Here the shaded bands represent the statistical uncertainty of the MadEvent integrations. 
(The statistical error is proportional to the cross section value and is therefore increasing 
for y/s > 2mt.) We see that Acrtf(A) is at the level of 1 — 2 fb for ^/s < 2mj and compatible 
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FIG. 10: Upper panels: the full relativistic cross section (t^''~''(A) obtained from MadEvent (red 
lines) and the LL (lower green dotted line), NLL (blue dash-dotted lines) and NNLL order (blue 
dashed lines) NRQCD cross sections. The uncut NNLL NRQCD cross section is also shown (upper 
green-dotted lines). Lower panels: difference of the relativistic and the NNLL NRQCD cross section 
(solid line), shaded bands represent the statistical uncertainty of the MadEvent integrations while 
the green dotted lines delimit the theoretical precision goal. We have chosen AMt = 15 for all 
panels on the left, and AMt = 35 GeV for those on the right. 



with zero within the statistical uncertainty of the MadEvent results for -y/i > 2mt. The 
greed dotted lines show the theoretical precision goal we have formulated in view of the 
experimental uncertainties expected at a future linear collider. The lines are determined 
from the quadratic sum of an energy-independent error of 5 fb and a 2% relative uncertainty 
with respect to the full NRQCD inclusive cross section prediction (including QCD effects) 
presented in Sec. |Vl Overall, we find that Aau^A) is much smaller than the theoretical 
precision goal indicated by the green dotted lines and compatible with the size of 0{v^) 
relativistic corrections. This demonstrates the excellent quality of the NRQCD phase space 
matching procedure we propose in this work and it should be adequate for the expected 
precision achievable at a future linear collider. 

In the upper panels of Fig. [TO]we have also shown the NNLL NRQCD cross section without 
the phase space matching contributions^ (upper green dotted lines). Except for the wave 



For = and without phase space constraints there are no NLL order corrections, and the LL and 
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function contribution in ^(C'y/yi which contributes less than 15 fb for the displayed 

energy range, this represents the approximation which has been used in all previous literature 
on the top pair threshold cross section where the top quark instability was accounted for by 
the replacement rule E ^ E + i Tf. We see that without phase space matching contributions 
the cross section is overestimated by about 30 — 50 fb at all energies. The relative discrepancy 
to the correct answer is particularly large for y/s < 2mt where the predictions that only 
account for the complex shift E ^ E + iTt do not vanish sufficiently fast. This unphysical 
feature is related to the fact that the top propagator of Eq. ([3]) leads to an overestimate of the 
contributions from NRQCD phase space regions with large top and antitop invariant masses. 
Thus the phase space matching corrections are essential to reach the intended precision goal 
for the theoretical predictions. 

Analysis for the full amplitude e~^e~ — )■ bbW^W~: 



We now consider the fully relativistic Standard Model amplitude for e~^e~ — ?■ bbW~^W~ at 
tree-level, a°'''^^{A). The results allow us to determine the remainder parts of the cross 
section, a^^^^{A). In the upper panel of Fig. [Tl]the fully relativistic cross section cr°''^°(A) 
obtained from MadEvent (red line) and the NNLL order NRQCD cross section (^Nmjcmi^ 
(blue dashed lines) are shown for AMf = 15 (left panels) and 35 GeV (right panels), where 
the parameter set of Eq. (|63l) has been employed and 10^ events have been used for inte- 
grations. The analytic formula for (^NRQCDti'"(^) J^^^ (EO]) with the electroweak top 
wave function contributions {iCyj^ ^^Zt,uni i^gjj^g replaced by the full interference coefficients 
^^v/A 1 froni Eqs. (122!) and (1231) . Up to relativistic corrections beyond NNLL order, which 
we neglect in the following discussion, the difference between the relativistic and the NNLL 
order NRQCD results is just a^^'^^^lA), which we use to determine the remainder contri- 
butions to the phase space matching coefficients Cy^^{A) and Cyj^{A). For the invariant 
mass cuts AMt = 15 and 35 GeV a^^^^{A) is displayed in fb units in the middle panel of 
Fig. [HJ The shaded region represents the statistical uncertainties of the MadEvent com- 
putation. For AMt = 35 GeV Crem^l^) ranges from to 5 fb and increases with the cm. 
energy. For AMt = 15 GeV it is almost energy-independent and below 1 fb for ^/s < 2mt 
and compatible with zero within the statistical MadEvent errors for -y/s > 2mt. The re- 
sults are fully consistent with a linear dependence in E, which conffims the structure and 
NRQCD counting of the {e^e~){e~^e~) forward scattering operators. Overall, the size of 
the remainder contributions is much smaller than the theoretical precision goal explained 
above and indicated by the green dotted lines. So the remainder part can be neglected 
for the theoretical predictions. With Eq. (126|) it is straightforward to fit for the remainder 
contributions to the phase space matching coefficients Cy^^ and Cyj^- We use an analysis 
for AMt = 15, 20, 25, 30 and 35 GeV, and assuming a hnear dependence on AMt the result 
of this fit reads 

AMt 



(7°(A) + C'°(A) 
C^^^''{A) + cTiA)) 



-1.2 ± 0.9) fb + (28 ± 6) fb 

rut 

-14 ± 25) fb + (260 ± 170) fb 

mt 



(64) 



NLL NRQCD cross sections agree. 
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FIG. 11: Upper panel: fully relativistic bbW^W~ cross section (t"''^'^(A) obtained from MadEvent 
with cut (red line) and without cut (green dotted lines), and NNLL order NRQCD cross section 
(T^^''(^^'^^^^ ( A) (blue dashed lines). Middle panel: af^^^{K) in fb units with the statistical uncer- 
tainties of the MadEvent computation (shaded region), and the theoretical precision goal (green 
dotted lines). Lower panel: A(Tnocut(A) = <7no™t ~ '^nrqcd (^) (black solid lines), and the sta- 
tistical uncertainties from MadEvent (shaded regions). We have chosen AM^ = 15 for all panels 
on the left, and AM^ = 35 GeV for those on the right. 



For the fit we added an energy- independent error of 1 fb to af^^'^{h) to account for higher 
order relativistic corrections oc {E/rrit)" with n > 1. 

Total cross section without cuts: 

Rather than predicting the inclusive cross section with cuts on the (anti)top invariant masses 
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it is also natural to ask about the total cross section for the case that no such cuts are applied 
on the bbW~^W~ final state. The phase space matching concept can be also applied in this 
case since the difference between the full theory phase space for the cases with and without 
invariant mass cuts constitutes again a hard NRQCD effect, much like the difference between 
applying cuts of different sizes. This can be easily understood because even without invariant 
mass cuts the phase space is constrained by kinematic bounds that in the tt threshold region 
(-y/i 2mt) also represent hard scales of order rrit. In the upper panels of Fig. [TT] the 
green dotted lines represent the full theory bbW~^W~ cross section without any cuts from 
MadEvent. The difference Acrnocut(A) = cr^J^^^ — cr^j^°cD^'"'"(A) is displayed by the black 
solid lines in the lowest panels of Fig. [TT]for AMt = 15 (left) and 35 GeV (right). In the 
figure the statistical uncertainties from MadEvent are again represented by the gray shaded 
regions. In contrast to predictions for the cross section with invariant mass cuts, here, A acts 
simply as a UV cutoff for the NRQCD phase space integration. Since the cross section is 
independent of the cutoff A, Ao'nocut(A) just compensates the A-dependence of the NRQCD 
cross section. Applying the same fitting procedure as for Eq. flM|) we obtain the following 
results for the remainder contributions to the phase space matching coefficients C*y/^ and 



'V/A 



C°(oo) + C'°(oo; 



rem 



(-18.9 ± 0.9) fb + (-42 ± 6)^^ fb 

mt 

(21 ± 26) fb + (90 ± 170) fb 



(65) 



Although we do not need the phase space matching coefficients Cyy^(oo) and Cyj^ (oo) for 
the predictions intended in this work, we give the results for future reference. 



IV. PHASE SPACE MATCHING WITH QCD EFFECTS 

In the previous sections we have discussed the concepts of the phase space matching 
procedure and computed the phase space matching coefficients neglecting QCD effects. We 
have demonstrated that the phase space matching coefficients can be reliably computed 
within NRQCD and that the remainder contributions which contain all effects that need 
the evaluation of full theory diagrams are small. For the precision goal we need to achieve, 
these remainder contributions can be neglected. We have also shown that for invariant mass 
cuts AMt = 15-35 GeV, power-counting breaking terms that arise from insertions of higher 
order relativistic operators do not spoil the nonrelativistic expansion. 

In this section we extend the phase space matching procedure to account for QCD cor- 
rections. Due to the additional QCD-induced interactions the required operator structure 
is more complicated. Up to NNLL order only the {e^e~){e~^e~) forward scattering opera- 
tors already discussed in the previous sections are required. At N^LL order we also need 
to account for imaginary phase space matching contributions to the Wilson coefficients of 
the {e~^e~){ti) top pair production operators to account for non-analytic energy-dependent 
terms in the matching relations. As for the previous sections we compute the phase space 
matching coefficients within NRQCD. Since for = we found the remainder contribu- 
tions to the coefficients to be negligible and since no kinematic or dynamical enhancement is 
expected for the QCD corrections, we ignore the QCD corrections to the remainder contri- 
butions in our analysis. For the determination of the QCD corrections we have to account 



35 



(a) A(A) 



'a(A) (&) 




(c) 'a(A) (d) ^(A) (e) A(A) 



FIG. 12: 0(q;s) QCD radiative corrections to e~^e~ — )■ tt — )■ bW~bW'^ originating from (a) time- 
like gluon exchange between 65-quarks, (6) time-like gluon exchange between tb and t6-quarks, (c) 
space-like gluon exchange between bb, and (d) interference between gluon radiation from b and b. 
Contribution (e) involving gluon radiation off the nonrelativistic top/antitop lines is suppressed 
with respect to the previous ones. The conjugated diagrams have to be added for each case. The 
definition of the top momenta pt in terms of the nonrelativistic momenta can be found in Fig. [2j 



for the QCD potentials and for the exchange of ultrasoft gluons among the top and antitop 
quarks and their decay products. Both types of corrections are analyzed in the following 
subsections. 



A. Ultrasoft Gluon Exchange and Nonperturbative Effects 

In the threshold region ultrasoft gluons, which carry momenta of order m^f^ ~ Ft, can 
propagate and lead to interactions among the top and the antitop quark and their decay 
products. These ultrasoft interactions are sometimes called "QCD interference effects" since 
the corresponding Feynman diagrams and their interference contributions all have compa- 
rable numerical size. The exchange of ultrasoft gluons starts to contribute at NLL order; 
the relevant diagrams at 0{as) in Coulomb gauge are displayed in Fi g. [121 I t is known that 



the QCD interference effects cancel in the total cross section at NLL |55l . |56| and NNLL or- 



der [27j if no cuts are imposed on the NRQCD phase space. Some of the interference effects 



have shown to cancel at NLL also for other inclusive observables where the top energy is 



integrated out, for example for the top quark three- momentum distribution [57[. For the top 
and antitop invariant masses Mt and Mj, on the other hand, the effects of ultrasoft gluon ex- 
change are essential because they directly affect the form of the distribution. The invariant 
mass distribution is further affected - at leading order - by nonperturbative effects. The pre- 
cise way how both effects enter the distribution depends on the details of the reconstruction 
prescription, see Refs. jH, H^] for a detailed treatment at large cm. energies ^/s ^ 2mt. The 
situation is, however, considerably simpler for the predictions of the inclusive cross section 
with the invariant mass cuts of Eq. (jl]) for AMt ^ F^. Because the inclusive cross section is 
an integral over Mt and in the resonance region, the main effects of the ultrasoft gluon 
exchange and of the nonperturbative contributions effectively correspond to a shift in AMt 
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of order Ft and Aqcd, respectively. This allows us to estimate that the size of the ultrasoft 
and the nonperturbative corrections to the phase space matching coefficients are of order 
asmtTtA"^ and mfAqcD/A^, respectively, with respect to the dominant NLL contributions 
of Eq. fl49l) . which are proportional to F^/A. The form of the nonperturbative corrections 
can also be understood from the point of view that there is an operator product expansion 
for the inclusive cross section. For the estimate for the ultrasoft gluons we have included a 
factor of as from the coupling of the gluon to the quarks. With Aqcd ~ c^s^t and A rrit 
both effects represent N'^LL corrections and are beyond the N^LL order level we consider 
in this work. Numerically they are suppressed by two orders of magnitude with respect 
to the NLL order phase space matching contributions. To confirm our argumentation for 
the ultrasoft corrections we explicitly compute in the appendix the phase space matching 
corrections to the {e~^e~){e'^e~) forward scattering operators that arise from the diagrams in 
Fig. O We use the invariant mass definition of Eq. (1331) and assume that the ultrasoft gluon 
can be resolved. While being not entirely realistic from the experimental point of view, this 
prescriptions should nevertheless give the typical size of the ultrasoft corrections. As shown 
in the appendix, we find that the ultrasoft phase space matching corrections indeed have 
the scaling behavior estimated above. In addition, they have a strong numerical suppres- 
sion that make them irrelevant. For the rest of this paper we therefore ignore ultrasoft and 
nonperturbative phase space matching corrections. 

B. Potential Interactions 

We now consider the phase space matching corrections originating from the potential 
interactions. Since the phase space matching procedure shares the properties of common 
loop graph matching computations, we find that once loop diagrams are involved, we also 
have to carry out phase space matching for subdiagrams in order to obtain local results for 
the matching coefficients. 

We start with the 0{as) contribution to the NRQCD cross section that arises from the 
insertion of one Coulomb potential. The diagrams to be considered are depicted on the LHS 
of the equality shown in Fig. [121 Following the approach of Sec. IIIIBl the result can be cast 
into the form 



iO{a) B 2 B 2 '^t^t f ^ "^i-E" " |(^1 + ^2) 

^RQCD =Nc ((C^r)' + i^ATf)^] dtldh + ^2r2)(^2 + ^2^2^ ^^(^1' ^2) + C.C. , 

A(A) 

(66) 

where in comparison to Eq. ( Hlj) we have to add the complex conjugate expressions because 
the computation involves a loop subgraph either to the left or to the right of the cut. For 
the Coulomb potential it is straightforward to derive 

A\t,,h) = v[Jm,E-^-l^^ , (67) 
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FIG. 13: 0{as) matching condition for the Wilson coefficients of e^e forward scattering opera- 
tors. 



where (g = |q|, v'f^ = - 
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(68) 



The term V is just the vertex diagram for ti production with one insertion of the Coulomb 
potential. Using Eq. fH6|) for the boundaries of the phase space integration we obtain the 
result 



Im[ln(— if) ] — 2 



A2 



8^2 m^r 



Re[iv] + O 



(69) 



0{as) Matching for the forward scattering operators: 

The first A-independent and non-analytic term in the brackets in Eq. (169|) is from the well 
known 0{as) contribution of the nonrelativistic Coulomb Green function, which contributes 
to the LL NRQCD cross section without cuts, see Eq. ([8]). The second term oc a^mtFi/A^ 
is local and contributes at NNLL order. Through the phase space matching procedure 
using Eq. the second term contributes to the imaginary Wilson coefficient Cy/A of the 
(e"'"e~)(e"'"e~) forward scattering operators Oy/A- The result reads 



C^/^(A)= 2zAr,(C^;Xi) 



2 

An . 



2a 



A2 



(70) 



and constitutes an 0{as) correction to the NLL matching contributions of Eq. fj49l) . The 
determination of the 0{as) contributions to Cy^^ arising from insertions of higher- dimension 
operators and potentials describing relativistic w^-corrections is carried out in analogy. Here 
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we just quote the final results and refer to Ref. 2^ for details: 



C'^;^(A)= 2ziV,(C^/X) 



2 mj 
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(71) 



where the definitions of the coefficients Cy°™j are given in Eqs. ( 120|) and ( l2Ti) . The coefficients 



Cl)), and Cy). come from one insertion of the potentials (k = p — p' 



'V/A 



'VIA 



V^^^P^ + P^^) 
1m\\? 



and 



mi 



(72) 



where we used V, 



-47ra and V, 



47ra/3, respectively. The origin of 



the other coefficients arises from insertions of the kinetic energy correction (kin), the time 
dilation correction (dil), the u^-suppressed S'-wave current 0^^2 (^^); the P-wave current 
C^p,3 (-P-wave) and the interference coefficient (int), as explained already after Eq. ( l53l) . The 
spin-dependent, momentum-independent potential vf^ |m?^ does not lead to a contribution 
at the order we consider. In analogy to our examinations for = we find that the phase 
space matching contributions originating from the relativistic insertions also contribute at 
NNLL order and are power-counting breaking. They also feature a relative factor h?/m1 
with respect to the terms in C'yy^(A) of Eq. (170]) . Through the numerical analysis of the 

C'yy^(A) coefficients which we carry out in Fig. [TU we find a similar situation as for the case 
of the corresponding tree-level matching coefficients C*^°^(A) studied in Sec. lIIIBl The phase 
space matching contributions in Cyy^(A) of Eq. ( ITOj) are about an order of magnitude larger 
than the ones from the relativistic corrections in Eqs. f l7T|) . which also cancel each other 



Since the potential Vs /m^ is momentum- independent, its insertion leads to a factorizcd expression for 
the correction to the Green fmiction that can only contribute to the phase space matching coefficient 5ci^ 
of the iS-wave tt current discussed in the following subsection. 
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FIG. 14: Absolute values of the contributions to the inclusive cross section from the 
0{as) phase space matching corrections Cyy^(A) of Eq. ([70l) (red lines) and Cy^^{A), i = 
-wave,int (red, blue, brown, magenta and cyan dashed lines, respectively), listed 
in Eq. (I7ip . The sum of the C'^y^(A) contributions is shown as the dotted black line. For compari- 
son we also display the tree-level NRQCD cross section without phase space matching contributions 
(black lines). The left panel corresponds to E = —5 while for the right E = 5 GeV. The values 
chosen for the input parameters can be found in Eq. ([63j) . and Ug = 0.1077. We have used the 
energy-independent form of the coefficients Cy°^- and Cyj^"^^ (i = 1,2,3) given by Eqs. (i20]l . 



'V/A,i 



partly due to their different signs. The cancellation is, however, less effective as for = 
due to the additional contribution coming from the Vr potential. Again, this does not apply 
to the contributions from the interference corrections in Cy^'J{A), which are comparable 

to Cyy^(A) for AMt > 25 GeV. As explained in Sec. IIIIBI before, this behavior is related 
to the large size of the interference coefficients Cyj^ ^ and not to power-counting breaking 
effects. As far as the size of higher order relativistic corrections are concerned we thus also 
find a good perturbative behavior for the 0{a.s) corrections to the phase space matching 
coefficients Cv/a{-^)- 

O(tts) Matching for the ti currents: 

Let us now discuss the third term in the brackets on the RHS of Eq. fl69l) which is oc Re[zi;]. 
It is non-analytic in the energy and can therefore not contribute to the phase space matching 
coefficients of the (e"''e~)(e"^e~) forward scattering operators OyJ^. This term contributes at 
N'^LL order and is also of theoretical interest as it illustrates how the phase space matching 
procedure is carried out at higher orders. We recall that the phase space matching follows 
the common rules of matching computations, where the matching involving diagrams with 
a higher number of loops first requires the matching of operators relevant for subdiagrams 
in order to deal with all terms that are non-analytic in the external momenta. To deal with 
the non-analytic term oc Re[u'] in Eq. dHH]) we need to determine the imaginary phase space 
matching term i6ci{A) contained in the Wilson coefficient Cy/^ i of the dominant 5- wave 
ti current operators Ov/A,p,i, see Eqs. (fTTl) . 

In the following we compute the contribution to the phase space matching term iSci{A) 
generated by an insertion of one Coulomb potential, as illustrated in Fig[T5l Suppressing the 
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FIG. 15: 0{as) matching condition for the Wilson coefficient of a tt production current. The 
black dot represents an insertion of a potential operator. 



external (anti)top quark spinors and the lepton current the amplitude for the cut diagram 
reads 



A(A) ' \2^P0 2mt ^ '' 2 J \2 ^ PO 2mt ^2 



^-.r*] (P-P' 



2mt 



(73) 



where ±p' is the 3-momentum of the top and the antitop quark in the c. m. frame, respec- 

(s) 

tively, and Vc = —AirCFas = — 4-7ra. For the determination of the matching condition for 
the ti current it is sufficient to consider on-shell external top quarks, and we therefore set 
p'^ = rritE, E > 0. By carrying out the angular integration and changing variables we 
obtain an intermediate result of the form (1441) with (|p| = {rritE — ^(ti + ^2))^''^) 



2n 



\p\y/mtE 



In 



VrritE + IpI 



iPl 



(74) 



The evaluation of the remaining integral with the invariant mass restrictions ( 1461) yields 

irrit ^ 1 + y/rritE 



t a 




E 



Re 



In 



rritV — \/rr^E 



8V2 mjTt „/ .m: 



+ 0{ 

Stt A3 V A5 



(75) 



The first term in the brackets is just twice the imaginary part of the vertex diagram I) 
defined in Eq. ( 168|) without phase space restrictions 



(76) 



and corresponds to the first term on the RHS of the matching relation in Fig. [151 The 
second term, on the other hand, gives a contribution to the matching term i5ci(A), and we 
obtain 



i5c[{k)= - (^^'■(A) - ^="(00)) = -ia 



Svr A3 ' 



(77) 



where the factor 1/2 is to compensate the factor 2 that appears in the optical theorem 
see Fig. [151 It is now straightforward to check that using the matching relation of Fig. 
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the phase space matching term i Sc\ (A) in the Wilson coefficient of the tt current operators 
Ov/A,p,i indeed accounts for the non-analytic term oc Re[2i;] in Eq. (169|) . This demonstrates 
the consistency of the phase space matching procedure at the loop level. We note that 
the same coefficient is obtained for the respective annihilation currents. Thus the imaginary 
phase space matching coefficients are not affected by the hermitian conjugation of the current 
operators in Eq. (ITBI) . 

The determination of the 0{as) contributions to i 5ci(A) arising from insertions of higher- 
dimension operators and potentials describing relativistic t;^-corrections is carried out in 



analogy. Here we just quote the final results and refer again to Ref. [29| for details: 



• r~kin,l/ a \ • ' 



i Sci (A) = —i a — , i 5c*' (A) = i a — , i Sc^™' (A) = —ia 



TT A ' ^ ^ ^ 37r A ' 1 ^ ^ 4^2 71 A 

t oci (A) = t a — , I oc^ [A) = I a — , t oc^ [A) = —z a — . (78j 

TT A Svr A vr A 

Note that the insertion of the current operators Ov/A,p,2 gives the mixing term i6cl con- 
tributing to the matching condition of Ov/A,p,i- An analogous mixing does not arise from 
insertions of the P-wave current Ov/A,p,3 due to angular momentum conservation. All terms 
shown in Eq. (1751) are of relative order A^/m^ with respect to the one for the dominant S- 
wave current in Eq. f l77|) . In the case of the interference correction this is also true because 



C. Perturbative Expansion and Higher Order Corrections 

In Sees. IIIIBI and IIVBI we have computed the and 0{as) phase space matching 

contributions to the {e~^e~){e~^e~) forward scattering operators. We have completed all 
contributions at NNLL order, and we have demonstrated that power- counting breaking 
terms do not spoil the nonrelativistic expansion. For the 0{as) contributions we have also 
shown that the Wilson coefficients of the {e^e~){ti) current operators receive phase space 
matching terms, which contribute at N^LL order. In this section we examine the a^-series of 
the phase space matching contributions. To render the analysis more transparent it is useful 
to distinguish the phase space matching contributions according to where they originate in 
the computation of the NRQCD cross section o"nrqcd(A). 

Contributions related to the imaginary part of the Coulomb Green function 

We ffist analyze the phase space matching contributions originating from the imaginary part 
of the Coulomb Green function Im[G'^] determined in Eqs. fH^ . fISU]) . f lTU]) and f l77|l . Since the 
imaginary part of the Coulomb Green function constitutes the leading order contribution of 
the factorization formula ( 126|) . these phase space matching terms represent the numerically 
dominant contributions. We define the NRQCD cross section associated to the imaginary 
part of the Coulomb Green function as 



A(A) 



(79) 
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where ti and ^2 are the invariant mass variables of Eq. (H2!) and fv,mt,u is the Coulomb vertex 
factor. The Coulomb vertex factor can be written in the form 



fv,mt,i'i\p\) 



nit 



G'M\P\). 



(80) 



Here G^^t v{^^ |p|) is the partially Fourier transformed Coulomb Green function with the 
first argument in position space at x = and the second in momentum space with the 
momentum p. Because of x = there is no dependence on the direction of p. At LL order, 
i.e. accounting only for iterations of the LL Coulomb potential of Eq. ([6]), ^(0, |p|) is 

known analytically: 



G^v,mtA^^ IpI) 



imt 1 
' Ak\^\ 1 - A 



k 

-.F,(2,l;2-A;-(l-ii^ 



(81) 



where k = ^y—mt{E + iTt) and A = Cpdsi'mti^) 'nit/'^k, and 2Fi{a^b; c; z) is the hypergeo- 
metric function. At NLL and NNLL order we use numerical results for the form factor fv,mt,v 
using the computational techniques we also emplo y fo r the determination of the Coulomb 



Green function at the respective orders, see Refs. |47l . |48| . These numerical methods have 



been applied and thoroughl y te sted before for predictions of the top three- momentum dis- 



tribution, see e.g. Refs. 
reduces to the form 



28|, |48[. Without phase space restrictions, i.e. for A = oo, Eq. ([79 



= 2N, 



+ (^^7)2) lm\G\a,v,mt,v) 



(82) 



In Fig. [16] we display Acr 



a 



NRQCD 



(A) -a. 



NRQCD 



[oo) for AMt = 15 GeV (left panel) 



and AMj = 35 GeV (right panel) using the LL (black solid line), NLL (red solid line) and 
NNLL (blue solid line) Coulomb Green functions, and the IS top mass scheme j28|, [Hsf with 
ntt = 172 GeV. The strong coupling is evaluated at the hard scale u = 1, as{mt) = 0.1077 
and the other input parameters are given in Eq. f[63]l . except for the electromagnetic coupling 
constant, which is also taken at the hard scale a^cdi^t) = 1/125.9. For simplicity we neglect 
here and in the following examinations of Sec. IIVCI the hard QCD corrections to the tt 
current matching coefficient, i.e. we set Ci(i^ = 1) = 1. The lines for Aa^"^^'^"^ represent 
the effects of the invariant mass cuts including iterations of the Coulomb potential to all 
orders. From our examinations in Sec. IIVBI we know that Ao"^™'*^"' exhibits hard and soft 
contributions that can only be separated by the phase space matching procedure. Thus 
unavoidably contains large logarithmic terms for any choice of the renormalization 
scale. The size of the higher order corrections discussed in the following should therefore 
be interpreted with some care as they might not reflect the quality of the Os-expansion 
when all logarithms are properly summed up. Despite this fact we see that the effects of 
including the 0{as) and 0{al) corrections to the Coulomb potential in Aa^^t*^"' show good 
convergence. For AMj = 15 GeV (35 GeV) the 0{al) corrections to the Coulomb potential 
(difference between the red and blue solid lines) lead to a shift between — 10 fb (—4 fb) and 
—5 fb (—2 fb). For AMt = 15 GeV this exceeds our theoretical precision aim (green dashed 
lines in the lower panel of Fig. [TO]) for y/s — 2mt < —2 GeV. The shift caused by the 0{a1) 
corrections to the Coulomb potential, however, quickly drops below the precision aim for 
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FIG. 16: Contributions to the inclusive cross section from the NLL, NNLL and N'^LL order phase 
space matching corrections in Acjp^jiP ^ (green, brown and black dashed lines, respectively) as a 
function of the total center of mass energy for AM^ = 15 (left) and AM^ = 35 GeV (right). 
The solid lines correspond to the phase space matching contributions Acr^™!*^''] obtained from the 
imaginary part of the exact cut Goulomb Green function at LL, NLL and NNLL (black, red and 
blue, respectively). 



AM^ > 15 GeV for all values of and we conclude that phase space matching corrections 
related to insertions of the 0{a1) corrections to the Coulomb potential, which contribute at 
N^LL order and higher, do not have to be considered. 

We now examine Aa'^™^'^''' computed from the phase space matching procedure carried 
out in the previous sections. Up to the N'^LL order it has the form 



Aa^"^!^^] =Im 

psm 



C'y{k)+C\{K) 

+ Im 



J NLL 



Im 



E 



(i),o 



J NNLL 



C^(A) + C1(A) + ^(c;.^^'^(A) + 4^)'°(A) 



i 5ci(A) G^{a, V, mt, v) 



(83) 



N3LL 



where we have indicated the NLL, NNLL and N^LL order corrections. The results for C?, 



^v/Ai ^vfl ^^"^ ^"^i have been given in Eqs. (l49l) . (TfOl) . ( l50l) . and ( j77l) . respectively. The 
N^LL corrections depend on the (MS renormalized) real part of the Coulomb Green function 
which at this order can be replaced by its LL expression see Eq. (jHl). Note that for actual 
NRQCD predictions the Green function G"^(a, mt, v) in Eq. ( l83l) has to be evaluated with 
~ to properly sum large logarithms. As mentioned above, we set z/ = 1 for the 
following examinations. In Fig. [16] the NLL and the NNLL approximations for Acr|iSiP ^ 
are displayed as the green and brown dashed lines, respectively. For AM^ = 35 GeV (right 
panel) we find that the NNLL corrections are about half the size of the NLL contributions. 
For AMj = 15_GeV (left panel) where we expect a worse Os-expansion due to the A- 
dependence of Cyjj^ ~ Ff/A and C^^^ ~ amtFj/A^, the NNLL corrections are only about 



15% smaller than the NLL contributions. The difference of Acp^*^ 



at NNLL order and 
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/\(jim[G=] accounting for the LL Coulomb potential (black solid lines) is between 3 and 7 fb 
for AMt = 35 GeV and between 7 and 25 fb for AMt = 15 GeV. For AMj = 15 GeV this 
exceeds our theoretical precision aim, visualized by the green dashed line in the lower panel 
of Fig. [TOl The difference is even larger with respect to Ao"^™'*^''' accounting for the NLL 
Coulomb potential (red solid lines). It is therefore required to also account for the full set 
of N^LL phase space matching contributions displayed in Eq. (!83|) . Unfortunately, at this 
time the full expressions for Cy^^ are unknown. They get contributions from two insertions 
of the leading Coulomb potential (contained in the LL Coulomb Green function) and from 
one insertion of the 0{as) correction to the Coulomb Green function (contained in the NLL 
Coulomb Green function). We have computed the contribution from two insertions of the 
Coulomb potential using the methods described in Sec. IIVBI The result reads 
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The logarithmic term is related to a NNLL order contribution to the imaginary anomalous 
dimension of Cy/A(A, z/). The fact that the logarithm vanishes for = A ~ reconfirms 
that the phase space matching contributions are hard effects. Although we do not have 
the complete result for Cy^^{A), using instead the result of Eq. (!84l) allows us to compare 

AcTp^f ^ at N3LL order (black dashed lines) with the numerical results for Aa^™^'-'"^ account- 
ing for the LL Coulomb Green function (black solid lines). We see that including the N3LL 
phase space matching contributions leads to a considerably improved agreement with the 
Coulomb resummed numerical results. For AMt = 15 GeV and 35 GeV the difference is 
always smaller than 5 fb, except when E < —5 GeV for AMt = 15 GeV. Since this is 
acceptable for our precision aim, we believe that the full set of N3LL phase space matching 
corrections should be adequate for the precision expected at a future linear collider. Since 
the full result for Cy^^{A) is unknown, we use for the time being as a substitute for the N3LL 

order terms in Adp^*^ ^ the numerical Coulomb-resummed expression of Eq. (179|) , evaluated 
with the NLL Coulomb vertex factor at the hard scale for as minus the NLL and NNLL 
order terms of Eq. (l83l) : 



psm 



Jn3ll =cr 



NRQCD 



(A, 1) - Im 



C'°(A) + C'°(A) 



— Im 



NLL 



J NNLL 



(85) 



Contributions related to the real part of the Coulomb Green function 

We now analyze the phase space matching contributions related to the real part of the 
Coulomb Green function. These phase space matching contributions are proportional to 
the type-1 imaginary Wilson coefficients iCyj^^ in Eqs. ( 122|) and ( 123|) . which describe the 

interference of the e'^e~ it ^ bbW^W~ diagram with diagrams for e~^e~ — )■ bbW^W~ 
with only either t or t at intermediate stages. As we have shown in Sec. IIIIBt the interference 
effects cause the largest phase space matching contributions among the (^(w^) relativistic 
corrections in the factorization formula f l26|) . and we therefore examine them separately. In 
analogy to the previous section we first define the Coulomb-resummed NRQCD cross section 
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FIG. 17: Contributions to the inclusive cross section from the NLL, NNLL and N^LL order phase 
space matching corrections in Acr^s'm' ^ (green, brown and black dashed lines, respectively) as a 
function of the total center of mass energy for AMt = 15 (left) and AMt = 35 GeV (right). The 
black solid lines correspond to the phase space matching contributions Aa^^^'-^''^ obtained from the 
real part of the exact cut Coulomb Green function at LL. 



with invariant mass restrictions arising from the real part of the Coulomb Green function: 
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Since the interference contributions are related to 0{v^) operator insertions in the factoriza- 
tion theorem ( 126|) . we only consider the vertex factor in the LL approximation as given in 
Eq. ( 180|) . The expression in the brackets involves the Coulomb-resummed generalization of 
the function A'°*'° given in Eqs. (135]) . The corresponding phase space matching contributions 
at N^LL order can be derived from Eq. f l26p and read 
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where we have again indicated the NLL, NNLL and N^LL order corrections. The results 

iint,0 /^int,l ^(l),int,0 i i:~int,l 

'V/A ■ 

respectively. The results for Cy^'^ are currently unknown. In the following analyses we 



for the C;™:,', C,VVr"'" and 6cf'^ have been given in Eqs. (El), (JTIl), (ESI) and (CHI), 



neglect them, i.e. we set them to zero. 

In Fig. [T7] the numerical LL Coulomb-resummed phase space effects Acr^®'*^'^' = 
(T^^qJj3(A) — (T^^qJj-)(oo) (solid black lines) and the corresponding NLL (green dashed lines), 
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NNLL (brown dashed lines) and N^LL (black dashed lines) approximations from the phase 
space matching procedure are shown. The left panel refers to AMt = 15 GeV and the right 
panel to AMj = 35 GeV. The strong and QED couplings are again evaluated at the hard 
scale ly = 1, i.e. as{mt) = 0.1077 and aqed('^t) = 1/125.9, and the other input parameters are 
given in Eq. (IMj) . As in the previous subsection we again neglect the hard QCD corrections 
to the tt current matching coefficient, i.e. we set Ci(z/ = 1) = 1. The situation we find is quite 
similar to the one discussed before for the imaginary part of the Coulomb Green function. 
The NLL order phase space matching contributions amount to around —13 (—20) fb, and 
the NNLL order contributions to about +18 (+15) fb for AMt = 15 (35) GeV. The NNLL 
corrections are quite sizeable and even exceed the NLL contributions for AMt = 15 GeV. 
It is conspicuous that the NLL and NNLL phase space matching corrections have opposite 
signs and cancel each other to a large extent. The sum of the NLL and NNLL phase space 
matching corrections differ from Aa^^^^"^ by at most 5 fb. The N^LL order phase space 
matching corrections, on the other hand, are much smaller than the NLL and NNLL order 
ones. Although we have neglected the Cy^^, which arise from diagrams with two insertions 
of the Coulomb potential, the difference of the phase space matching contributions up to 
N^LL order and the exact LL Coulomb-resummed result amounts to less than 2 fb for all 
AMt between 15 and 35 GeV. We conclude that the Cy^^ are numerically small and that 

the N^LL phase space matching contributions are more than adequate for our theoretical 
precision aim. 



Contributions related to the other 0{v'^) relativistic corrections 

As the third class of phase space matching contributions we now examine the corrections 
that arise from insertions of 0{v'^) suppressed operators other than the interference cor- 
rections just discussed above. Up to N'^LL order the contributions of these matching cor- 
rections to the inclusive NRQCD cross section have the form {i = {kin, dil, t;^, P-wave}, 
j = {r, s, kin, dil, v"^, P-wave}, k = {r, s, kin, dil, w^}) 
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where we have again indicated the NLL, NNLL and N^LL order corrections. The results for 
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^vm'° ^^"^ ^^i^ have been given in Eqs. fl5^ . flTTjl . fl55|l and fl78|) . respec- 
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tively. As for the interference coefficients the results for the Cifjj^ are currently unknown 
and are neglected in the following analyses. 



In Fig. Ha AcTpsl^ ^ is displayed at NLL ( green dashed lines), NNLL (brown dashed line) 
and N^LL (black dashed lines). It is striking that the NLL contributions are at the sub-fb 
level and an order of magnitude smaller than the NNLL terms. The small size of the NLL 
terms is, however, due to an almost complete cancellation in the sum of the individual C'y/^ 
coefficients (see Fig. [S]). The size of the individual coefficients is at the level of 1 — 3 fb 
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FIG. 18: Contributions to the inclusive cross section from the sum of the phase space matching 
arising from insertions of 0{v^) relativistic corrections other than the interference contributions 
at NLL (green), NNLL (brown) and N^LL (black) order, as a function of the total center of mass 
energy for AMt = 15 (left) and AMj = 35 GeV (right). 



and comparable in size to the individual NNLL coefficients Cy/^^ (see Fig. [Ml) . The NNLL 
coefficients, on the other hand, do not cancel to the same extent due to the additional 
contribution from the potential Vr'^\ see Eq (I7T]) . Thus the situation concerning the a^- 
expansion for the third class of the phase space matching contributions is similar to the 
other classes of phase space matching contributions discussed above. This is also confirmed 
by the size of the N'^LL corrections which amount to at most 1 to 2 fb. Although we do not 
have a numerical Coulomb-resummed expression for a more thorough comparison, the results 
indicate that the a^-expansion is also well under control for the third class of phase space 
matching contributions. Due to the overall small size of these contributions we conclude 
again that keeping the phase space matching contributions up to N'^LL order is more than 
adequate to reach our theoretical precision goal. 



V. NUMERICAL ANALYSIS 

In Sec. nil CI we have demonstrated for the case ag = that the phase space matching 
contributions are essential in order to compensate for the fact that the previous NRQCD 
predictions for top threshold production overestimate the full theory cross section by a sub- 
stantial amount. The problem of the previous NRQCD predictions is related to the fact 
that upon shifting the energy hj E ^ E + iTt, in order to account for the top quark fi- 
nite lifetime, the NRQCD tt phase space becomes unrestricted and extends to unphysical 
regions as a consequence of the nonrelativistic expansion. The phase space matching proce- 
dure removes these unphysical phase space contributions and implements the information on 
possible experimental cuts into the NRQCD predictions. These phase space matching con- 
tributions are represented by imaginary contributions to the Wilson coefficients of NRQCD, 
see Sec. [TTl 

In this section we analyze the complete set of N'^LL phase space matching contributions 
determined in the previous sections for predictions of the inclusive NRQCD top pair thresh- 
old cross section with cuts on the invariant masses of the top and the antitop quark defined 
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in Eq. (IMI) . We also compare the size of the phase space matching contributions to the other 
and previously known types of electroweak effects relevant for the top pair threshold cross 
section: the NNLL QED corrections, the hard electroweak corrections 22| and the (type-1) 



finite lifetime corrections, which are not related to phase space restrictions [27|. Together 



with the phase space matching contributions, which we call type-2 finite lifetime corrections, 
these three classes of contributions constitute all effects of the electroweak interactions on 
the threshold cross section. 

To start we collect all phase space matching contributions to the inclusive NRQCD cross 
section up to N^LL order. In the previous sections we have for simplicity neglected the hard 
QCD and QED matching corrections contained in the Wilson coefficient ci{u = 1) of the 
leading {e~^e~){tt) top pair production operator. Accounting for these matching corrections 
the complete set of N'^LL phase space matching contributions can be derived from the 
factorization formula in Eq. (12^ and takes the form 
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where we have specifically indicated by brackets the NNLL and N^LL order contributions. 
Explicit expressions for the NLL and NNLL hard QCD / QED matching conditions h^i^ and 
hf'^ are given in Eq. ( !T9|) . The terms A(j^'''"^(A, 1) are the N*^LL phase space matching 
contributions to the inclusive cross section with the h^^''^^ set to zero. The result for Aa^^^ 
reads (i = int, dil, kin, f ^, P-wave) 
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where expressions for Cy^^ and the Cy^^ have been given in Eqs. (149|) and ( 154|) . respectively. 
The term Aa'^^'"'" has the form [j = r, s, int, dil, kin, v"^, P-wave) 
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where expressions for Cy/j^ and the C*;^/^ have been given in Eqs. (fTOj) and ( ITT]) , respectively. 
Finally, Ao"^^^^ reads {i = {int, dil, kin, w^, P-wave}, j = {r, s, int, dil, kin, f^, P-wave}, /c = 



We do not discuss here the effects of the e+e luminosity spectrum since it is determined for the most 
part from experimental measurements and simulations. 
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where {Aa^Ja ')nJ'll is a numerical expression defined in Eq. fl85|) . and the results for Cyj^''^ 
and i^c'^'^ have been given in Eqs. and (frS]) . respectively. The results for the C*;^/^ are 
currently unknown. For the analysis we carry out in the following we set them to zero. 
Together with an analytic determination of {Aa^^ ^)n3ll we plan to compute them in a 

separate publication. 

In the following analysis we use 



172 GeV for the top mass in the IS mass scheme [28|, 
and all matching coefficients are evaluated at the scale rrit {v = 1). For the QCD 
coupling we use as{mt) = 0.1077 and for the QED (MS) coupling aqcdifnt) = 1/125.9. All 
soft matrix element contributions are evaluated for the QCD and QED couplings at the 
velocity renormalization parameter u = 0.2 which corresponds to yUgoft = 34.4 GeV for the 
soft and to /iusoft = 6.88 GeV for the ultrasoft scales. For the evaluation of the hard one- 
loop electroweak corrections we choose mniggs = 130 GeV. All other parameters are given in 
Eqs. dSS]). 

In Fig. [19] we show the phase space matching corrections to the inclusive ti threshold 
cross section at NLL (black dotted lines), NNLL (black dashed lines) and N^LL order (black 
solid lines), the NNLL QED corrections (blue lines), the NNLL hard one-loop electroweak 
corrections (green lines) and the type-1 finite lifetime corrections (red lines) as a function of 
the cm. energy ^/s. The left panel shows the results for an invariant mass AMt = 15 GeV 
and the right panel for AMt = 35 GeV. The QED effects arise from the electromagnetic 
correction to the QCD Coulomb potential (see text after Eq. (I7j)) and the one-loop QED 
matching correction to the Wilson coefficient ci of the tt current, see Eq. f lTTl) . The hard one- 



loop electroweak corrections are encoded in the coefficients CyJ^^ also shown in Eq. ([T 



The result for CyJ^^ has been obtained in Ref. 122 



V/A, 

The type-1 finite lifetime corrections 
represent all finite lifetime corrections which are not related to phase space constraints. 
They consist of the corrections generated by the imaginary interference matching coefficient 
iCyj^ 1' see Eq. ( ITTl) . the time dilation corrections to the Green function shown in Eq. ( I30ll 
and the contributions from the renormalization group summation of phase space logarithms 
contained in the coefficients Cy/A of the {e~^e~){e^e~) forward scattering operators given in 
Eq. ( l3T|l . The matching and time dilation corrections are known at NNLL order and the 
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FIG. 19: Sizes of the different contributions to the inclusive cross section arising from electroweak 
interactions as a function of the total center of mass energy for AMj = 15 GeV (left) and AMt = 
35 GeV (right): (green line) NNLL hard one- loop electroweak effects from Ref. [2^, (red line) 
NNLL finite lifetime corrections from Ref. [271 ] . (blue line) NNLL QED effects, and phase space 



matching corrections at NLL, NNLL and N'^LL (dotted, dashed and solid black lines, respectively). 
The blue dotted lines correspond to the expected experimental uncertainties at the LC. 



summation of phase space logarithms at NLL order j27j. The QED, hard electroweak and 
type-1 finite lifetime corrections do not depend on phase space restrictions and are therefore 
identical in both panels. In Fig. [19] the blue dotted lines represent a rough (and likely 
optimistic) estimation of the expected experimental uncertainties at a future linear collider 
consisting of an energy-independent error of 5 fb and a 2% relative uncertainty with respect 
to the full prediction, both being added quadratically. 

We see that the QED (blue lines) and the type-1 finite lifetime corrections (red lines) are 
sizeable (at the level of 40 fb) only in the peak region just below ^/s = 2mt. Above and 
below the peak region the QED corrections are quite small and do not exceed 5 fb. Above 
and below the peak the type-1 finite lifetime corrections amount to —15 to —10 fb. Due to 
their different signs the QED corrections and the type-1 finite lifetime corrections cancel each 
other to a large extent in the peak region. The hard electroweak corrections (green lines) 
represent a multiplicative factor of -1.2% to the total cross section and are therefore very 
small below the peak and at the level of 12-13 fb above the peak region. We see that the 
phase space matching contributions represent the largest of the four classes of electroweak 
effects. In contrast to the other classes of electroweak effects they do not decrease strongly 
for energies below the peak region. For AMt = 15 GeV the N^LL phase space matching 
contributions amount between —85 and —65 fb and for AMt = 35 GeV they are between 
—45 and —35 fb. The overall size of the phase space matching corrections decreases for 
larger values of the top invariant mass cut AMt. We emphasize, however, that the results 
obtained in this work are valid only for moderate values of AMt in the region between 15 and 
35 GeV. For invariant mass cuts below 15 GeV the phase space constraints are not related 
anymore to hard effects and for invariant mass cuts substantially above 35 GeV matching 
contributions that need to be computed from full theory diagrams have to be included. The 



The small size of the hard electroweak eorrections displayed in Fig. [12] is obtained for the QED coupling 
defined at the scale of the top mass mt. 
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FIG. 20: Total inclusive top pair production cross section from NRQCD: starting from the pure 
QCD NNLL prediction (black dashed line), we add step- by-step the QED corrections (blue line), 
the hard electroweak corrections (green line), the type-1 finite lifetime corrections (red line) and 
the N'^LL phase space corrections (black solid line) for AMt = 15 GeV. 



relatively flat behavior of the phase space matching contributions is related to the fact that 
the dominant phase space matching contributions are energy-independent. The small linear 
dependence on ^/s is related to the ^/s dependence of the virtual 7 and Z propagators of 
the basic e+e~ — j- tt process and the peak-like structure comes from an imaginary phase 
space matching contribution to the {e~^e~){ti) top pair production operator which enters 
the N^LL inclusive cross section in terms of a time-ordered product. This time-ordered 
product leads to a non-analytic dependence on the energy, see Eqs. fl26|l and fl92|) . In Fig. [19] 
we have also displayed the phase space matching contribution to the inclusive cross section 
at NLL (black dotted lines), NNLL (black dashed lines) and N^LL order (black solid lines) 
in order to show the convergence of the phase space matching procedure. The results show 
that the expansion related to the phase space matching procedure is particularly good for 
larger values of AMt and still well under control for AMt = 15 GeV. We note that the 
rather small size of the NNLL corrections (difference of black dotted and dashed lines) for 
AMt = 15 GeV arises from a cancellation between different independent NNLL corrections, 
see Sec. HVCl 

In Figs. |20] and [21] the size of the four different types of electroweak corrections is shown 
for predictions of the total inclusive cross section. Starting from the pure QCD NNLL 
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FIG. 21: Total inclusive top pair production cross section from NRQCD: starting from the pure 
QCD NNLL prediction (black dashed line), we add step- by-step the QED corrections (blue line), 
the hard electroweak corrections (green line), the type-1 finite lifetime corrections (red line) and 
the N'^LL phase space corrections (black solid line) for AMt = 35 GeV. 



cross section (black dashed lines), which accounts only for the electroweak effects from 
the basic e~^e~ Z ^ ti process and the LL finite lifetime effects through the energy 

replacement rule E ^ E + iTt), we add step-by-step the QED corrections (blue lines), the 
hard electroweak corrections (green lines), the type-1 finite lifetime corrections (red lines) 
and the N^LL phase space corrections (black solid lines). In Fig. [20] we have AMt = 15 GeV 
and in Fig. [21] we have AMt = 35 GeV. Only the phase space matching corrections depend 
on AMt. We again see that the phase space matching contributions exceed by far the 
other electroweak corrections. Since we have already discussed the size of the individual 
types of electroweak corrections in our analysis of Fig. [191 we concentrate here only on the 
behavior of the predictions for energies below the peak region where the cross section is 
small. Here the phase space matching corrections are very large and amount to changes of 
more than 50% percent for AMt = 15. These large corrections are related to the unphysical 
phase space contributions contained in the pure QCD prediction which are a consequence 
of the nonrelativistic expansion and the energy replacement rule E ^ E + iTt- We stress 
that these unphysical effects cannot be cured by adding more of the higher order QCD 
corrections because they originate from modifications to the nonrelativistic tt phase space 
products caused by the top width. Thus in order to obtain realistic predictions for the top 
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threshold cross section it is essential to account for the phase space matching corrections. 



VI. CONCLUSION 



The effects of the finite top quark lifetime are an essential ingredient for predictions 
of the top pair production rate close to threshold y/s ~ 2mt, where rrit is the top quark 
mass. Because the top width Fj has approximately the same size as the typical top quark 
kinetic energies, finite lifetime effects already enter the leading-order predictions and cannot 
be treated as corrections. An important consequence is that the top pair cross section is 
only well-defined with a given set of prescriptions how the observed top decay final states 
are accounted for in the cross section. This entails that the cross section also accounts for 
non-tt processes which lead to the same observed final state, and that the cross section 
can depend on experimental cuts on kinematic variables such as the reconstructed invariant 
masses. When matching to the NRQCD effective theory these finite lifetime effects can 
be integrated out for observables that are inclusive on the top and antitop decays, and 
lead to imaginary matching contributions to the Wilson coefficients of NRQCD operators. 
These imaginary matching coefficients are much like the complex indices of refraction in the 
Maxwell theory of light propagation in an absorptive medium. The cross section including 
the finite lifetime effects can then be obtained from the absorptive part of the e"'"e~ — )■ e"'"e~ 
forward scattering amplitude using the optical theorem. 

One can distinguish two types of imaginary NRQCD matching coefficients. The type-1 
contributions 27| account for the (Cutkosky) cuts through the top decay final states in full 
theory diagrams. They lead e.g. to the well known quark bilinear top width term and also 
account for the interference effects mentioned above. Insertions of the associated opera- 
tors also cause UV divergences in the nonrelativistic ti phase space integrations since the 
resulting unstable top propagator i/(po ~ P^/2mt + iTt/2) lifts the stable particle disper- 
sion relation po = P^/2?7it and allows for arbitrarily large final state top invariant masses. 
These UV divergences require the introduction of {e~^e~){e~^e~) forward scattering opera- 
tors, which acquire an imaginary anomalous dimension and sum large logarithms of the top 
velocity in the tt final state phase space. The type-2 imaginary matching contributions are 
the matching conditions of the Wilson coefficients of these {e~^e~){e'^e~) forward scattering 
operators. They encode the information on the experimental cuts used for the cross sec- 
tion measurement. Since the type-2 matching contributions are not related to hard virtual 
fluctuations, but to real final state configurations with large top quark off-shellness, we call 
them phase space matching contributions. 

In this work we have determined and analyzed the type-2 phase space matching condi- 
tions up to N^LL order in the nonrelativistic expansion for cuts on the reconstructed top and 
antitop invariant masses Mf, of the form |M( f — mt| < AM^, with AM^ between 15 and 
35 GeV and neglecting the W boson width. We have demonstrated that the numerically 
dominant effect of the phase space matching is to remove the unphysical NRQCD phase 
space contributions that do not pass the cut. This is because the nonrelativistic unstable 
top propagator i/{po — p^/2mt + iTt/2) overestimates by far top and antitop fluctuations 
with large off-shellness. The remaining numerical contributions to the phase space matching 
conditions coming from calculations of relativistic full theory diagrams were found to be well 
below 5 fb for the cross section, which is negligible for the experimental precision one can 
expect at a future linear collider. From the field theoretic point of view, the procedure of car- 
rying out the phase space matching agrees with the common matching and renormalization 
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methods for stable particle theories. Thus at higher orders in the nonrelativistic expansion 
it is required to account for the phase space matching contributions of subdiagrams to re- 
move non-analytic matrix element terms from the matching equations and to achieve that 
the matching coefficients are analytic in the external energy. In the phase space matching 
procedure at N'^LL order for top pair production at threshold one has to also include the 
phase space matching for (e+e~)(tf) top production operators. 

Since the phase space matching procedure we have carried out involves the computation 
of NRQCD phase space integrations with a hard cutoff related to AMj, our results contain 
power-counting breaking contributions. This means that phase space matching contributions 
coming from insertions of higher order operators can give contributions that are formally 
lower order. We have shown that such power-counting breaking contributions do not spoil 
the nonrelativistic expansion and that power-counting breaking can be ignored from the 
practical point of view. As far as the expansion is concerned, we have found that the 
N^LL (0(a^)) corrections to the phase space matching contributions need to be determined 
to meet the experimental precision expected at a future linear collider. At this time these 
N^LL corrections are not yet fully known analytically, and their determination shall be 
addressed in subsequent work. 

Our final numerical results have been given in Sec. |Vl Leaving aside the effects from 
the e^e~ luminosity spectrum, which are known to distort the cross section shape and 
normalization in a quite substantial way, we find that the phase space matching contributions 
to the cross section exceed by far the other types of electroweak corrections, which are known 
from previous work. The phase space matching contributions are between —85 and —35 fb 
for invariant mass cuts AMf between 15 and 35 GeV and are essential for realistic theoretical 
predictions. In the peak and the continuum region {y/s > 2mt) they amount to 6 to 10%. 
They are particularly important in the region below the peak {^/s < 2mt) where the cross 
section decreases and the unphysical off-shell contributions of the NRQCD ti phase space 
become dominant. Here the phase space matching contributions can amount to more than 
50%, and they ensure that the cross section has the correct physical behavior. 

Phase space matching is also important for predictions of the top pair threshold cross 
section if no kinematic cuts are imposed, since the NRQCD phase space contributions in- 
volving off-shell top and antitop quarks still lead to large unphysical contributions. Here the 
phase space matching contributions from relativistic full theory diagrams are numerically 
important and cannot be neglected (Sec. IIII CD . In this work we have determined these full 
theory contributions for = 0, i.e. at NLL order. We finally note that the phase space 
matching procedure can also be carried out using exclusively full theory computations with 
kinematic cuts to determine the imaginary type-2 matching coefficients. This approach does 
not involve any power-counting breaking contributions and allows to determine the phase 
space matching contributions more easily for AMf > 35 GeV. At this time these full the- 
ory computations are only known for = 0, which allows to carry out this phase space 
matching approach at NLL order. To go beyond the NLL level the results for the 0{as) 
and 0{al) corrections to the e~^e^ cross section are required for the final states that arise 
in top pair production. Such results are not available at this time. 

Acknowledgment s 

This work was supported in part by the DFG Sonderforschungsbereich/Transregio 9 
"Computergestiitzte Theoretische Teilchenphysik" and the EU network contract MRTN- 



55 



CT-2006-035482 (FLAVIAnet). P.R. thanks M. Beneke for useful discussions. Feynman 
diagrams have been drawn with the packages Axodraw 59 | and Jaxodraw 6oI . 



Appendix A: QCD interference effects 

In this appendix we compute the 0{as) ultrasoft corrections to the imaginary phase space 
matching coefficients Cy/A of the {e~^e~){e~^e~) forward scattering operators Oy/A for the 
invariant mass prescription explained in Sec. IIV A[ The corresponding diagrams in Coulomb 
gauge are shown in Fig. [T^l For the interference diagram [T2i we need to define the top and 
antitop invariant masses in the presence of an additional ultrasoft gluon in the final state. 
As a toy prescription that can be easily implemented analytically we assume that we can 
resolve the gluon down to an infrared scale A. For gluon energies larger than A we define 
the top and antitop 4-momentum as the sum of 4-momenta of their decay products, bW^ 
and bW~, respectively. The top and antitop invariant masses are then defined exactly as in 
Eq. ( p3|) . The result for the ultrasoft phase space matching corrections obtained from this 
prescription should also be generic for the typical size of corrections for other, more realistic 
invariant mass prescriptions. We can write the 0{as) contributions to the imaginary phase 
space matching coefficients Cy/A from each of the diagrams in Fig. [T^-d as {i = a, b, c, d) 

^V/A,(i) = '^^c{Cv'/A,lY ^(7) . (Al) 

The resulting contributions to the inclusive cross section read 

AaJ^^ = iV. [{C^jn' + iC^rf] c'l!^ ■ (A2) 

The computation of the ultrasoft corrections for energies close to threshold is performed 
using NRQCD Feynman rules for the top/antitop propagators and ultrasoft gluon couplings. 
The gluon momentum k is neglected when appropriate according to the scaling k ~ rritv'^. 
In Coulomb gauge the time-like gluon propagator has the form z/k^ and the transverse 
propagator is —kV/k'^)/{k'^ + ie). We cut the diagrams as indicated by the red dashed 
lines in Fig. [12] using the well-known Cutkosky rules for the transverse gluon propagator 
and Eq. fl39|l for the top and antitop propagators. We note that our computation is quite 
similar to the one presented some time ago in Ref. [56|, where, however, no phase space cuts 
were considered. The main conclusion in the work of Ref. 5™ was that for the total cross 



section (i.e. without phase space restrictions) the contributions from the ultrasoft diagrams 
in Fig. [12] cancel. This serves as an important cross check of the computations we carry out 
here. 

The contribution from the time-like gluon exchange between the bb pair (Fig. [T2k plus 
the conjugated diagram) reads: 

us _ 4.„ 4 /• d^P f d'k (3|A(fc°,k)p-|B(fc°,k)p) 

^ ]. h c c (A3) 

(ti - imtTt){t2 - imtTt){ti - 2mtk° + imtTt){t2 + 2mtk^ + irritTt) 



We disagree with the corresponding resuU given in Ref. [56| with respect to the sign of the term | B (fc" , k) ^ . 
However, this does not ahcr the conchision that the contribution cancels in the absence of phase space 
restrictions. 
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where the invariant mass variables ti, t2, are defined in Eq. (H2|) . The same variables are 
employed for the computation of the diagrams in Figs. [T2b and c. In Eq. (1A3P we have 
already performed the integrations over the bW~^ and bW~ phase space variables, which 
yields the functions 



B(A;°,k) = 2 



gs^t 



In 





k 






k 


-%e) 



gsTt 2x 



Ikl 2x + l 



2Ty 



In 





k 


— ie 


+ 


k 


— ie 



k 



(A4) 



where x = M^/m'f. It is straightforward to check that Eq. ( ]A3I) vanishes if no bounds are 
imposed on the integration over the top energy: carrying out the integration by residues 
one obtains a purely imaginary number, which cancels out when adding the conjugate di- 
agram. This confirms that the diagram vanishes if there are no cuts on the phase space 
integration. For the invariant mass cuts of Eq. (H6l) we proceed by performing the integra- 
tion over the 3-momentum and the energy of the virtual gluon. After the trivial integration 
over the p angles we obtain a representation of the time-like gluon exchange diagram of the 
form 



27r3 



dtidt2 



A(A) 



^JmtE-{h + t2)/2 



^2) • 



(A5) 



We use this generic form for all the QCD interference diagrams, i = a,b, c, d. Note that 
expression ( lASP is compatible with Eq. (144|) . i.e. the A(ti,t2) functions in both expressions 
have the same normalization. For the time-like gluon exchange from diagram [T2b we obtain 



^ \ (2x-\ 



9 V2X + I7 ) (tl + ^2)2 + 4(mtrj)2 



: I (2^1^)' + t2(ti + t2)) arctan + ^(ti - t^) In 



{tl ^t2}\. 



(A6) 



The time-like gluon exchange between th and t6. Fig. WDo . is computed analogously: 



16 i Cpgsrri 



(27r)4 ] (27r) 



1 



(tl - 2mtk^ + imtTt){t2 + 2mtA;0 + zmjFt 



tl ^ t2 
A;0 ^ -A;0 



c.c. 



(A7) 



The terms shown explicitly in Eq. (]A7ll represent the contribution from the tb gluon exchange 
(first diagram in Fig. [T2b). The contribution from the ib gluon exchange is obtained with 
the replacements ti <H- t2, and — )■ —k^, as indicated in Eq. flA7p . For this contribution 
we find it more convenient to perform first the A;°-integral. Again, it is easy to check that 
if the ^^-integration is done by residues, one obtains a result which is purely imaginary, so 
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this contribution vanishes for the total cross section without phase space restrictions. The 
expression can be cast into the form of Eq. (1A5|) with 



1 /2x - 1 



9 V 2x + 1 



2\ -1 



^(a)(il, ^2) 



(A8) 



Let us now turn to the space-hke gluon exchange between the final state bb pair in Fig. [T2h . 
The result reads 



where 



(ti - 2mtk^ + imtVt){t2 + 2mtk'^ + imtVt) 



c.c. 



(A9) 



C(A;°,k) 



fif^Ft 2x - 1 



Ikl 2x+l 



2k2 



In 



k'>- 


k 


— ie 




k 


— ie 



(AlO) 



The expression can be cast into the form of Eq. (lASP with 



A(c)(tl, ^2) 



Ci?as /2a; — 1 



In 



tl + (mfEt 

127r \2x+lJ [V \ (2mt)4 
(ti-t2)(4mjri + f (ti+ta))/ 



An h r ^ 

— arctan h |ti ^ t2| 

3 mtEt 



(ti-t2)(ti+t2-fmir 

^ (ti + t2)2 + 4Kri)^ 



^2 

arctan arctan 

V mtVt rritTt 



t) I'tl + {nitVt 
m 



2 In 



{2mty 



(All) 



The |k| integration in c^^-^ yields an infrared divergence, which we have regularized with the 
cutoff A mentioned at the beginning of this appendix. This IR divergence is cancelled by 
a corresponding IR divergence in the real gluon emission diagram in Fig. [T2l i. as shown 
below. Attaching a space-like ultrasoft gluon to the top or antitop lines such as in diagram 
Wlk yields an additional v factor, so the corresponding contributions are suppressed in the 
nonrelativistic expansion with respect to those in Figs. [T^-d. This way we only need to 
consider real gluon emission from the bottom quark lines. Due to the additional gluon in 
the final state the relation between the variables ti_2 and the nonrelativistic loop momenta 
(p", p^) has to be modified. For the momentum routing displayed in Fig. [T2ji, ti and ^2 are 
defined as 



V2 ^ 2m, / 



2m+ 



E 



P 



2mt 
2mt 



(A12) 
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The result from diagram 
4 



-id) 



then reads 
d^p f d^k C(|k|,k)2 



1 



1 



(ti + 2mt|k| + imtTt){t2 + 2mt|k| — irritT 



{ti - imtVt){t2 + imtVt) 
+ c.c. . 



(A13) 



For this contribution the |k|-integration extends to values such that the phase space factor 
IpI = ^JmtE — mt|k| — {ti + t2)/2 remains a real number. As anticipated above, the |k|- 
integration is IR-divergent and we introduce the cutoff A as in the case of diagram [T2b. The 
expression can be cast into the form of Eq. flASP with 



CpOis f 2x — 1 



4 In 



+ 



127r V2x+ 1, 
2{t2-imtTt) h{t2) 



2mtE - (ti + t2] 



mi 



2 In 



_ti -t2 + 2imtTt h{ti + t2 



In 



h{t2) + h{ti+t2) 
h{t2) -h{ti+t2) 



(2mi)2 
{h ^ t2} 



c.c. 



(A14) 



where 



h{y) = ^/mtE - {y - imtTt)/2. 
h{y) = ^mtE-y/2. 



(A15) 



^us 

-(d)- 



We explicitly see the cancellation of the infrared divergent terms in the sum of c"^^ and c" 
Again it is easy to check that after integrating over the top energy without restrictions, 
both contributions cancel completely. 

Expanding the ultrasoft phase space matching contributions nitE^ TritTt <^ A^ up to terms 
of order 1/A^, we obtain 



^/2m 



16) 



'(c) 



37r2 



Cpas 1 



1 /2x 



9\2x 



1 3 1^ + 1^ 2 - 1 - - sinh"Vr 



37r^ 



V2', 



m; 



37r3 



A3 



2x - 1 
2x + l 



A2 



A3 [ A2 
ln2- 1 - V2-sinh"Vr 



A 



Amt 
2A2 



72^ 



-id) 



37r3 



2x-iyE,r ^E, 

2x + l A 1 A 



rritE 
3A^ 



97r 



1 ^tEf 
3^ A2 



(^V2 + sinh"^(l)jln 



Am* 



rri^k2_ 
A2 97r 



2A2 
(A16) 
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with 



ki = rfi + 5v^-6V21n2-^-i^ + 7r(31n2-i) + 7r2(v^+^-21n2) 
+ (5 _ - 91n2) sinh-i(l) + ^ [sinh-i(l)]' - SLi^l^^—^) , 

k2 = d2 + ^{5-12\n2) + ^. (A17) 

The constants di and d2 have been evaluated numerically and read 

di = -4.961 , ^2 = -17.75. (A18) 

Adding up the results of all diagrams we are - as anticipated from the general arguments 
discussed in Sec. IIV Al - left with a correction to the cross section of order mtF^/A^ with 
respect to the NLL phase space correction: 



i=a,b,cA 



A- = ^f^V^(vr^fl + ^-ln2 + sinh-^(l))-A:,-A:a. (A20) 



where 



277r2 V2x + ly A2 

Apart from the cancellation of the In A infrared divergences already pointed out above, we 
also find that all logarithms of Fj cancel in the sum of all diagrams. This is expected be- 
cause the phase space matching contributions represent hard effects. A proper evaluation of 
the ultrasoft phase space matching corrections therefore also requires to employ the strong 
coupling in Eq. ( IMUj) at the hard scale. For nit = 172 GeV, Mw = 80.425 GeV and 



Cpasimt) = 0.1436 we find that A"*^ = 0.004 rrifFf/A^. We thus find that the ultrasoft 
corrections based on our invariant mass prescription have an additional strong numerical 
suppression factor. Although this result might not be generalized to other invariant mass 
definitions, it nevertheless supports the conclusion that ultrasoft effects are in general irrel- 
evant at the level of precision we aim for in this work. 
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